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ABSTRACT ■  A  preliminary  analysis  of  the  dispersion  process  of  liquid 
droplets  into  an  aerosol  cloud  is  presented.  This  assessment  of  the 
complex  phenomena  is  not  all  encompassing,  but  is  believed  to  shed  some 
light  upon  the  behavior  of  the  various  mechanisms  which  govern  the  ex¬ 
pansion  of  explosively  dispersed  liquids.  This  analysis  does  provide 
a  basic  understanding,  and  provides  data  for  preliminary  design  esti¬ 
mates  . 

The  equations  presented  are  based  on  the  concept  of  a  fluid  piston 
in  the  first  three  of  four  regimes  of  motion.  This  fluid  piston  is 
essentially  a  gas  bubble  surrounded  by  a  liquid  shell.  As  the  gas 
bubble  expands  it  forces  the  liquid  outward  to  form  the  characteristic 
explosive  dispersion  cloud.  The  fluid  piston  degenerates  into  discrete 
droplets,  whose  motion  is  governed  by  aerodynamic  forces  in  Regime  IV. 

These  equations  of  motion  were  developed  from  experimental  data 
collected  with  two  different  size  devices.  The  agreement  with  this 
data  indicates  that  scaling  requirements  are  satisfied. 
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FOREWORD 


This  investigation  was  conducted  during  the  period  May  1967  until 
December  1968  in  direct  support  of  development  programs.  This  effort 
was  funded  under  AirTask  A35-532-050/216-1/W11-62X  and  additional  fund¬ 
ing  was  provided  under  Marine  Corps  Commitment  Authorization  WR-27-8- 
7721. 


This  report  summarizes  the  first  part  of  a  proposed  two  part  effort 
-=to  develop  a  understanding  of  the  mechanics  of  explosive  dispersion  of 
liquids.  The  Part  1  effort  consists  of  a  dispersion  model,  describing 
the  mechanisms  by  which  the  liquid  is  explosively  forced  into  an  expand¬ 
ing  cloud. 
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INTRODUCTION 


The  dispersion  of  an  aerosol  cloud  of  liquid  droplets  by  means  of 
the  detonation  of  an  explosive  charge  within  a  canister  of  liquid  is 
becoming  of  increasing  interest  for  military  purposes.  Furthermore, 
the  process  may  have  potential  c'vi.’ian  applications  for  the  dispersion 
of  liquids  (or  even  solid  powders'  for  such  purposes  as  firefighting 
or  fertilizing  crops,  in  which  ch:^icals  might  be  dispersed  rapidly 
over  large  areas. 
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In  order  to  properly  employ  explosive  dispersion  techniques  with 
maximum  effectiveness,  we  must  better  understand  the  complex  physical 
phenomena  making  up  the  explosive  dispersion  process.  In  particular, 
we  have  to  understand  those  factors  controlling  the  expansion  of  the 
cloud  in  the  atmosphere  so  that  the  devices  can  be  designed  to  disperse 
the  material  into  the  desired  cloud  configuration  with  the  desired  con¬ 
centration  . 
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In  this  report  we  have  attempted  a  preliminary  analysis  of  the  dis¬ 
persion  process  based  on  the  experimental  information  which  is  available 
and  which  can  be  published  in  an  unclassified  report.  We  don’t  in  any 
sense  represent  this  as  "the"  analysis  of  the  cloud  dispersion  process, 
but  do  believe  it  will  shed  some  light  on  the  complex  physical  phenomena 
which  govern  the  expansion  of  explosively  dispersed  clouds  and  enable  us 
to  at  least  begin  to  predict  the  behavior  of  these  clouds.  Considerably 
more  experimental  and  analytical  work  will  be  required  before  we  can 
claim  fully  to  understand  the  process. 
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BACKGROUND 

A  typical  device  for  the  explosive  dispersion  of  a  liquid  is  illus¬ 
trated  in  Fig.  1.  It  consists  of  a  metal  canister  into  which  the  ends 
are  welded.  Usually  the  canister  is  scored  with  longitudinal  grooves 
placed  approximately  15  degrees  apart  to  insure  that  breakup  of  the 
canister  is  uniform.  A  metal  tube  which  contains  the  explosive  burster 
charge  is  located  along  the  axis  of  the  canister  and  is  constructed  so 
that  the  explosive  charge  may  be  inserted  after  the  canister  is  filled 
with  liquids  and  sealed.  A  detonator  is  placed  in  close  proximity  to 
the  explosive  charge  as  shewn.  Table  1  lists  the  dimensions  and  capac¬ 
ities  of  tv’o  typical  devices ,  one  containing  three  pounds  of  liquid 
and  the  other  83  pounds. 

Upon  the  detonation  of  the  burster  explosive  the  liquid  is  forced 
outward  with  such  energy  that  it  breaks  up  into  aerosol  droplets  which 
disperse  in  a  few  milliseconds  to  form  a  toroid-shaped  cloud  as  shown 
in  Fig.  2.  The  concentration  of  liquid  droplets  is  much  less  near  the 


FIG.  1.  Typical  Three-Pound  Explosive  Dispersion  Device. 
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TABLE  1.  Dimensions  and  Capacities  of  Explosive 
Dispersion  Devices. 


Three-pound 

Eighty-three  pound 

Liquid  weight 

3  lb 

83  lb 

Liquid  density  (typical) 

1.80  slug/ft3 

1.80  slug/ft3 

Liquid  volume 

5.18  x  10-2  ft3 

1.43  ft3 

Explosive  weight 

0.053  lb 

0.728  lb 

Explosive  density 

3.10  slug/ft3 

3.10  slug/ft3 

Explosive  volume 

5.70  x  10"4  ft3 

7.29  x  10~3  ft3 

Canister  length 

0.833  ft 

1.60  ft 

Canister  diameter 

0.333  ft 

1.12  ft 

Canister  vail  thickness 

0.065  in 

0.062  in 

Canister  material 

Alum.  (6061  T-4) 

Steel  (1020) 

Burster  tube  volume 

11.4  x  10-4  ft3 

1 

8.11  x  10~3  ft3 

Burster  tube  diameter 

0.625  in 

1.065  in 

Burster  tube  wall 
thickness 

0.035  ir. 

0.035  in 

Burster  tube  material 

Steel  (4130) 

Steel  (4130) 

End  plate  thickness 

0.25  in 

0.375  in 
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center  of  the  cloud  than  In  the  portion  near  the  outer  edge.  The  shape 
of  the  cloud  is  complex  in  that  "spikes"  of  liquid  are  forced  outward 
from  the  main  body  of  the  cloud  as  the  expansion  proceeds. 

In  previous  investigations  it  has  been  assumed  that  the  cloud  is  a 
toroidal  mass  in  which  the  radial  velocity  is  constant  with  respect  to 
circumferential  position.  Attempts  were  made  to  plot  the  edge  of  the 
cloud  as  it  expanded  radially  and  it  was  implicitly  assumed  that  the 
edge  of  the  cloud,  as  seen  from  the  side,  was  composed  of  the  same  drop¬ 
lets  throughout  the  expansion  process.  Robbins,  for  example,  in  Ref.  1 
and  2  photographed  the  outward  expansion  of  the  cloud  and  obtained  a 
composite  radius  versus  time  plot  for  representative  devices.  The  pre¬ 
sent  authors,  along  with  Douglas  Kinney,  used  the  same  technique  in  an 
early  report  and  then  attempted  to  derive  an  empirical  relation  for  the 
equation  of  motion  based  on  the  Gurney  formula  for  the  expansion  of  ex¬ 
plosive  products  (Ref.  3) . 

It  soon  became  apparent  in  the  present  study,  however,  that  such  an 
approach  was  not  going  to  yield  meaningful  results  in  terms  of  under¬ 
standing  the  physical  phenomena  which  were  controlling  the  expansion  of 
the  droplets  of  liquid.  The  cloud  does  not  expand  uniformly  with  res¬ 
pect  to  circumferential  position.  It  is  well  known  that  spikes  of  liquid 
follow  the  fragments  or  strips  of  canister  material  as  they  travel  out¬ 
ward.  Those  "strip  spikes"  travel  outward  considerably  further  than 
the  main  body  of  the  cloud.  Furthermore,  another  type  of  spike,  which 
we  call  a  "squirt  spike",  occurs  at  the  time  of  canister  rupture.  The 
squirt  spikes  are  caused  by  liquid  squirting  outward  through  the  cracks 
in  die  canister  as  it  ruptures  along  the  longitudinal  grooves.  The 
squirt  spikes  initially  travel,  some  200  ft/sec  faster  than  the  cloud 
body,  but  slow  down  more  quickly. 

An  observer  viewing  the  cloud  from  the  side  will  see  first  the 
squirt  spikes  leading  the  remainder  of  the  cloud.  Later  the  cloud  body 
will  pacc  tne  squirt  spikes  and  finally  the  droplets  composing  the  strip 
spikes  will  move  out  in  front.  Thus,  the  expansion  sequence  is  complex 
in  that  different  physical  processes  govern  the  expansion  of  different 
portions  of  the  cloud.  It  certainly  may  be  possible  to  develop  a  single 
empirical  equation  which  can  be  used  to  describe  the  position  of  the 
cloud  at  any  time,  but  such  an  equation  would  not  assist  us  in  under¬ 
standing  the  physics  of  the  dispersion  process.  To  do  this  we  must 
treat  each  portion  of  the  cloud  separately  and  attempt  to  develop  equa¬ 
tions  of  motion  which  described  the  physical  processes  which  are  occur¬ 
ring  in  the  cloud. 


EXPERIMENTAL  DATA 

Only  a  limited  amount  of  detailed  experimental  data  of  a  quantita¬ 
tive  nature  is  available,  and  most  of  this  is  for  the  experimental  three- 
pound  device.  We  have,  therefore,  based  this  study  primarily  on  the 
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information  available  for  the  three-pound  device.  We  have  also  attempted 
to  scale  up  the  analysis  to  match  information  which  is  available  for  a 
device  containing  83  pounds  of  liquid. 

Inspection  of  photographs  taken  at  approximately  6000  frames  per 
second  by  means  of  a  Fas  tax  camera  showed  that  during  the  dispersion 
of  a  cloud  from  a  three-pound  canister  a  shock  wave  breaks  away  from  the 
canister  at  about  the  time  the  walls  begin  to  move  outward.  This  shock 
wave  has  been  observed  by  other  experimenters.  It  is  assumed  that  the 
shock  wave  is  the  result  of  the  detonation  wave  from  the  burster  charge 
-propagating  outward  through  the  liquid  and  the  canister  walls  and  into 
the  atmosphere.  The  shock  wave  is  quite  weak  and  has  an  initial  speed* 
of  only  200  or  300  ft/sec  greater  than  the  speed  of  the  canister  walls 
as  they  expand.  Thus  the  majority  of  the  energy  of  the  detonation  wave 
has  been  given  up  to  the  liquid  and  the  canister  walls  before  the  wave 
emerges . 

Examination  of  photos  of  the  same  devices  taken  with  a  Cordin  camera 
at  framing  rates  of  approximately  125,000  frames  per  second  indicates 
that  the  canister  expands  before  breaking.  The  time  from  first  motion 
of  the  canister  until  rupture  is  about  0.0004  seconds.  The  canister 
begins  to  break  near  the  center  of  the  longitudinal  grooves  and  liquid 
immediately  begins  to  squirt  out  through  the  resulting  orifices  which 
are  formed.  This  phenomenon  can  be  seen  in  Fig.  3,  which  is  a  reproduc¬ 
tion  of  a  typical  Cordin  camera  framing  sequence.  The  liquid  which 
squirts  outward  from  the  canister  as  it  is  breaking  into  strips  is  that 
which  forms  the  squirt  spikes. 

As  the  canister  expands  further  it  breaks  into  discrete  fragments 
or  strips  which  fly  outward.  The  lid  and  the  bottom  of  the  canister 
are  also  separated  and  fly  up  and  down  respectively,  although  at  lower 
speeds . 


The  main  body  of  the  cloud  can  begin  to  expand  into 
only  after  the  canister  has  broken.  Since  the  fragments 
behind  the  squirt  spikes,  and  are  traveling  more  slowly, 
are  not  visible  in  photographs  taken  of  the  side  of  the 


the  atmosphere 
are  initially 
they  initially 
cloud.  The  for¬ 
ward  edge  of  the  cloud  is  the  forward  edge  of  the  squirt  spikes  for  the 

photographs  of 
itrates  the  ex- 
The  framing 


first  few  milliseconds.  Figure  4  is  taken  from  overhead 
a  typical  three-pound  device  during  dispersion  and  illusl 
pansion  pattern  of  the  spikes  and  main  body  of  the  cloud, 
rate  was  approximately  3400  frames  per  second. 


The  fragments  do  not  decelerate  as  rapidly  as  the  droplets  com¬ 
prising  the  squirt  spikes  and,  as  the  cloud  expands,  the  fragments 
catch  up  with  and  pass  the  squirt  spikes.  Because  the  fragments  create 


bpeed  is  used  throughout,  as  opposed  to  velocity,  since  the  dir¬ 
ection  is  only  assumed  to  be  radial  and  speed  represents  k.ie  gross 
magnitude  of  all  the  particles  collectively. 
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wakes  behind  them  they  drag  a  certain  portion  of  the  liquid  droplets 
behind  them  as  they  fly  outward.  Those  droplets  form  the  strip  or 
fragment  spikes,  and  those  droplets  pass  the  squirt  spikes  and  become 
the  le.ading  edge  of  the  cloud  as  it  expands. 

The  droplets  comprising  the  main  body  of  the  cloud  expand  at  a 
slower  rate  than  the  droplets  in  either  of  the  spikes.  We  estimate, 
however,  that  the  droplets  in  the  main  body  comprise  more  than  90%  of 
the  total  mass  of  the  cloud.  Hence  it  is  the  rate  and  distan:e  of  ex¬ 
pansion  of  the  main  body  droplets,  not  the  spikes,  which  is  important. 
Fortunately  the  spikes  are  thin  and  tenuous  and  in  the  side  view  the 
main  body  of  the  cloud  can  usually  be  differentiated  from  the  spikes. 
Thus,  the  information  on  the  expansion  rate  and  distance  of  the  main 
body  of  the  cloud  usually  can  be  obtained  from  a  careful  assessment  of 
the  experimental  data.  The  main  body  of  the  cloud  never  expands  to  the 
extent  that  the  fragment  spikes  do,  but  we  consider  that  the  actual  use¬ 
ful  diameter  of  the  cloud  is  that  of  the  main  body,  and  not  that  of  the 
spikes  . 


SHOCK  WAVE  PROPAGATION  DISTANCE  AND  RATE 

The  shock  wave  from  a  three-pound  device  was  photographed  using 
Fas  tax  cameras  at  framing  rates  of  approximately  6000  frames  per  second. 
Figure  5  illustrates  the  shock  wave  radius  as  a  function  of  time  for  a 
typical  three-pound  device.  The  initial  shock  wave  speed  appears  to  be 
about  1500  ft/sec  as  recorded  by  both  cameras.  The  shock  decays  slowly, 
but  must  eventually  reach  a  minimum  speed  of  approximately  1140  ft/sec, 
the  speed  of  sound  in  air.  It  was  impossible  to  read  the  shock  wave 
position  beyond  approximately  10  milliseconds  because  it  became  too  weak 
to  be  visible  against  the  backdrop. 

The  only  quantitative  shock  speed  information  we  have  available  is 
that  for  the  three-pound  device.  The  information  gained  from  different 
experimental  firings  agrees  quite  well.  It  has  not  been  feasible  to 
attempt  to  obtain  shock  speed  information  for  larger  devices  since  the 
detonation  of  the  larger  bursters  would  damage  the  experimental  gridded 
backdrop  which  is  used  to  read  the  shock  position. 


FRAGMENT  RADIUS  AND  SPEED 


The  intial 
obtained  from  th 
shown  in  Fig.  6. 
with  the  detonat 
by  a  tenth  of  a 
approximately  th 
mately  1100  ft/s 
initial  speed  of 


radius  versus  time  relation  of  the  fragments  could  be 
e  Cordin  camera  photographs  for  two  typical  firings  as 
Because  of  the  difficulty  of  synchronising  the  cameras 
ions,  the  two  sets  of  data  are  displaced  from  each  other 
millisecond.  Nevertheless  the  slopes  of  the  curves  are 
e  same  and  indicate  initial  fragment  speeds  of  approxi- 
ec.  The  initial  speed  of  the  fragments  is  also  the 
the  main  body  of  the  cloud  since  the  liquid  (other 
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than  that  In  the  squirt  spikes)  follows  the  fragments  outward.  Hence, 
for  the  three-pound  device  at  least,  the  initial  speed  of  the  main  body 
of  the  cloud  is  about  1100  ft/sec.  The  fragments  appear  to  decelerate 
slowly  as  indicated  by  Fig.  7,  which  gives  typical  fragment  and  cloud 
distances  as  a  function  of  time  for  the  first  0.010  seconds.  These  data 
were  taken  from  Fas  tax  camera  photographs  at  a  framing  rate  of  about 
6000  frames  per  second.  Note  that  the  initial  average  fragment  speed 
must  be  approximately  1100  ft/sec  (as  indicated  by  the  Cordin  data)  for 
the  fragments  to  reach  the  point  of  emergence  at  the  point  and  time  in¬ 
dicated  on  the  Fastax  data.  This  occurs  at  approximately  3.5  feet  and 
3  milliseconds.  Information  is  not  available  on  the  initial  fragment 
speeds  for  the  larger  devices . 


SQUIRT  SPIKE  RADIUS  AND  SPEED 

The  squirt  spikes  have  an  initial  speed  greater  than  that  of  the 
fragments  since  they  are  projected  outward  through  the  cracks  in  the 
canister  because  of  a  pressure  differential  between  the  inside  and  out¬ 
side  of  the  canister.  Experimental  data  plotted  from  Cordin  camera 
frames  is  illustrated  in  Fig.  8.  It  can  be  seen  that  for  the  three- 
pound  device  at  least,  the  initial  speed  of  the  squirt  spikes  is  about 
1225  ft/sec,  or  some  100  to  150  ft/sec  faster  than  for  the  fragments 
of  the  main  body  of  the  cloud. 

It  is  hard  to  differentiate  between  the  squirt  spikes  and  the  main 
body  of  the  cloud  during  the  first  several  milliseconds  of  flight.  Ac¬ 
cording  to  the  Fastax  data  shown  in  Fig.  9  and  10,  which  are  intended 
to  be  for  the  main  body  of  the  cloud,  the  initial  speed  of  the  cloud  is 
about  1250  ft/sec.  This  compares  well  with  the  Cordin  data  for  the 
squirt  spikes.  This  is  because  the  main  body  cannot  be  read  during  the 
first  few  milliseconds  and  we  see  only  the  squirt  spikes  on  the  Fastax 
frames . 


MAIN  BODY  RADIUS  AND  SPEED 

Out  tc  about  3  milliseconds  and  3.5  feet  radius  the  position  of  the 
main  body  can't  be  determined  precisely  and  must  be  assumed  to  lie  some¬ 
where  between  the  fragments  and  the  squirt  spikes.  Some  of  the  main 
body  is  undoubtedly  trapped  behind  the  fragments  but  some  of  the  main 
body  droplets  probably  are  forced  outward  between  the  fragments  because 
of  the  internal  pressure  of  the  expanding  detonation  products.  Thus 
the  main  body  will  have  a  speed  i  .  'his  period  of  between  1100  and 
1300  ft/sec.  Beyond  three  milliseconds  the  main  body  can  be  differ¬ 
entiated  quite  well  and  its  position  can  be  plotted  as  shown  in  Fig.  9 
and  10.  Because  the  difference  in  speed  between  the  fragments  and  the 
squirt  spikes  is  not  great  (at  least  for  the  three-pound  device)  we  have 
assumed  the  initial  speed  of  the  main  body  to  be  that  of  the  squirt  spikes, 
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FIG.  7.  Fragment  and  Cloud  Radius  Three-Pound  Device 
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Experimental  data  for  speeds  as  a  function  of  volume  for  the  main 
body  of  the  cloud  are  shown  in  Fig.  11.  These  data  were  obtained  by 
graphically  measuring  the  slope  of  the  tangent  to  the  faired  radius 
versus  time  curve  of  Fig.  9  at  various  points.  The  speed  was  plotted 
as  a  function  of  cloud  volume  in  order  to  facilitate  comparison  with 
the  analytical  equation  of  motion,  which  was  more  readily  derived  as  a 
function  of  cloud  volume.  The  volume  is  assumed  to  be  equal  to  the  cloud 
area  times  the  canister  length.  For  the  purposes  of  the  analysis  the 
cloud  was  assumed  to  remain  constant  in  thickness  during  its  expansion. 

We  know  this  is  not  true,  but  the  error  introduced  by  the  assumption  is 
not  great  and  the  assumption  simplifies  the  analysis. 

The  main  body  radius  versus  time  curve  for  a  typical  83-pound  device 
is  shown  in  fig.  12.  It  can  be  seen  that  the  initial  cloud  speed  appears 
to  be  much  greater  than  for  the  three-pound  device.  This,  we  believe,  is 
a  result  of  a  much  greater  initial  squirt  spike  speed  for  the  83-pound 
device.  The  initial  slope  of  the  curve  appears  to  be  on  the  order  of 
2000  ft/sec.  This  is  not  unreasonable  if  we  assume  that  the  steel  canister 
of  the  83-pound  device  might  rupture  before  expanding  as  far  as  the  aluminum 
canister  of  the  three-pound  device.  If  this  were  the  case  the  pressure 
differential  between  the  inside  and  outside  of  the  steel  canister  would  be 
much  greater  than  for  the  aluminum  canister  at  the  time  of  rupture.  Hence 
the  squirt  speed  would  be  greater  for  the  steel  canister. 


As  will  be  shown  later  an  initial  main  body  speed  of  approximately 
1380  ft/sec^  is  more  reasonable  and  when  that  value  is  used  the  two  de- 
vi.  compare  quite  well  when  allowances  are  made  for  scaling.  Figure  13 
shows  the  main  body  speed  versus  cloud  volume  data  for  the  83-pound  de¬ 
vice  in  which  an  initial  main  body  speed  of  1380  ft/sec  is  assumed. 


ANALYSIS  OF  CLOUD  DISPERSION 


PHYSICAL  MODEL  OF  CLOUD  DISPERSION 

In  order  to  derive  equations  of  motion  for  the  expansion  of  the 
explosively  dispersed  cloud  it  was  necessary  to  develop  a  simple  model 
of  the  dispersion  process.  We  assumed  that  the  device  was  composed  of 
a  cylindrical  metal  canister  filled  with  liquid  as  shown  in  Fig.  1.  A 
cylindrical  burster  tube  extends  through  the  center  of  the  canister  and 
is  partially  filled  with  explosive.  The  presence  of  the  metal  of  the 
burster  tube  is  neglected,  and  the  liquid  is  assumed  to  fill  the  re¬ 
mainder  of  the  canister  with  no  allowance  for  ullage. 


This  is  the  initial  speed  as  calculated  in  a  later  section. 
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CLOUD  VOLUME  (FT)® 


FIG.  11.  Cloud  Speed  as  a  Function  of  Volume 
Three-Pound  Device. 


CALCULATED 

+  EXPERIMENTAL  DATA 


VOLUME  (FT*) 


FIG.  13.  Cloud  Speed  as  a  Function  of  Volume  83  Pound  Device 
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The  burster  explosive  is  assumed  to  detonate  instantaneously  and 
expand  to  fill  the  remainder  of  the  burster  column  with  detonation  pro¬ 
ducts.  The  detonation  products  then  expand  radially  outward.  A  deton¬ 
ation  wave,  preceding  the  detonation  products,  is  assumed  to  move  radi¬ 
ally  outward  from  the  burster  tube,  becoming  a  shock  wave  as  it  enters 
the  liquid.  The  cylindrical  shock  wave  propagates  radially  outward  from 
the  burster  column  until  it  ente-s  the  atmosphere  surrounding  the  can¬ 
ister.  But,  since  the  strength  of  the  shock  wave  emerging  from  the 
canister  is  small,  we  assume  that  the  energy  of  the  detonation  is  given 
up  entirely  to  the  liquid  and  fragments. 

The  liquid  is  assumed  to  expand  radially  outward.  As  noted  pre¬ 
viously,  we  assume  that  the  cloud  does  not  increase  in  thickness  during 
the  expansion  process,  but  remains  equal  in  thickness  to  the  length  of 
the  canister. 

To  simplify  the  analysis  we  have  divided  the  dispersion  process  into 
several  regimes,  each  of  which  will  be  treated  separately.  This  is  nec¬ 
essary  because  the  physical  processes  governing  each  regime  are  different 
so  that  it  is  only  reasonable  to  expect  that  different  equations  of 
motion  will  be  required.  We  have  attempted  to  develop  equations  of 
motion  which  will  grossly  predict  the  behavior  of  the  cloud  without  undue 
complexity.  As  we  gain  further  understanding  of  the  phenomena  involved 
it  will  be  possible  to  refine  or  rework  the  analysis  to  provide  a  more 
exact  model. 


INITIAL  REGIME 

The  Initial  Regime  is  assumed  to  begin  when  the  burster  explosive 
detonr.tes.  During  this  regime  the  burster  explosive  is  transformed  from 
solid  explosive  filling  a  portion  of  the  burster  tube  to  detonation  pro¬ 
ducts  at  very  high  temperature  and  pressure  completely  filling  the  bur¬ 
ster  tube. 


Knowing  the  detonation  pressure  and  the  initial  volume,  which  is 
assumed  to  be  the  initial  volume  of  the  solid  explosive,  we  can  deter¬ 
mine  the  pressure  of  the  detonation  products  at  the  end  of  the  Initial 
Regime  by  means  of  the  poly  tropic  equation  for  the  gas: 


P  «  P  — 

og  dg  l  V 


where 


n  =  polytropic  exponent 

P  »  gas  pressure  at  the  end  of  the  Initial  Regime 
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P,  -  detonation  pressure  of  the  explosive 
dg 

V  =  gas  volume  at  the  end  of  the  Initial  Regime 
°g 

V,  “  solid  explosive  volume  at  detonation 
dg 

According  to  Ref.  4  the  polytropic  exponent  is  about  3.0  at  the  time  of 
detonation  hut  decreases  rapidly.  We  will  assume  an  average  value  of 
n  ™  2.0  for  the  Initial  Regime.  The  detonation  pressure  of  a  typical 
explosive,  from  Ref.  4,  is  87,000  atm  or  1.84  x  10  psf. 

Three-Pound  Device 

For  a  typical  three-pound  device  the  solid  explosive,  as  can  be 

.seen  from  Table  1,  fills  half  the  burster  tube  so  that  V  =  2V,  .  Sub- 

og  dg 

stituting  in  Eq.  1  we  get 

P  -  (1.84  x  108)  (0.5)2  ■=  4.6  x  IQ7  psf 
°g 

Eighty-Three  Pound  Device 

The  detonation  pressure  and  the  average  value  of  n  were  assumed 
th  same  for  the  83-pound  device  as  for  the  three-pound  device.  From 
Ta-le  1  the  remaining  initial  conditions  are: 

V,  =  7.29  x  10"3  ft3 
dg 

V  -  8.11  x  10~3  ft3 
og 

Substitution  of  the  above  values  in  Eq.  1  results  in  a  final  pressure: 

P  *  1.49  x  10  psf 

Aft  r 


This  pressure  is  three  times  that  of  the  three-pound  device  because  the 
free  volume  in  the  burster  tube  for  the  83-pound  device  is  proportionately 
much  less  than  in  the  three-pound  device. 


REGIME  I 

f 

*  Regime  I  begins  when  the  burster  tube  begins  to  expand,  forcing  the 

l  I  liquid  and  canister  walls  outward.  We  assume  that  the  shock  wave  from 

j  |  the  detonation  imparts  negligible  speed  to  the  liquid,  which  is  forced 

f  outward  solely  by  the  pressure  of  the  expanding  detonation  products, 

f  I  The  canister  walls  are  forced  outward  in  turn  by  the  liquid.  In  this 


regime  the  gaseous  detonation  products  are  forced  to  expand  against  the 
inertia  of  the  liquid  and  the  canister  walls,  the  resistance  of  the 
stretching  canister  walls,  and  the  pressure  of  the  surrounding  atmosphere. 

We  assume  that  the  canister  walls  expand  so  that  the  canister  remains'1 
cylindrical  and  its  length  remains  constant.  It  is  further  assumed  that 
no  liquid  vaporizes  during  the  process. 

When  the  metal  of  the  canister  reaches  its  ultimate  elongation,  the 
canister  ruptures  along  the  longitudinal  grooves,  signifying  the  end  of 
Regime  I. 

In  this  regime  we  assume  our  system  (or  control  volume)  to  be  the 
liquid  cylinder  which  is  being  dispersed.  It  is  bounded  on  the  inside 
by  the  expanding  detonation  products  and  on  the  outside  by  the  walls  of 
the  cylinder  and  the  ends  of  the  canister  (see  Fig.  14). 


Equation  of  Motion 


The  energy  equation  for  a  non-flow  system  is  the  basis  for  the  equa¬ 
tion  of  motion.  We  assume  that  no  mass  is  transferred  across  the  system 
boundaries  and  that  heat  transfer  and  changes  in  potential,  internal, 
and  chemical  energy  are  negligible.  The  only  work  done  on  or  by  the  sys¬ 
tem  is  expansion  work  so  that  the  energy  equation  becomes 


-  6W 


exp. 


dE. 


(2) 


where  is  the  average  kinetic  energy  of  the  liquid  mass.  It  is  assumed 
to  be  given  by  the  relation 


dE,  **  m.v  dv 

k  la 


(3) 


where 


m^  ■  mass  of  the  liquid  constant 

v  =  average  speed  through  the  liquid  shell 
a 

The  expansion  work  is  composed  of  three  parts;  the  work  done  on  the 
liquid  by  the  expanding  detonation  products,  the  work  done  by  the  liquid 
in  expanding  the  metal  canister,  and  the  work  done  by  the  liquid  in  push¬ 
ing  aside  the  atmosphere.  In  this  regime,  since  the  work  done  in  expand¬ 
ing  the  canister  is  much  greater  than  the  atmospheric  work,  the  latter  cn.n 
be  neglected. 

The  work  done  on  the  system  by  the  detonation  products  is  given  by 
the  relation 
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W  =  -  IP  dV 
g  J  g  g 


where 


=  pressure  of  the  gaseous  detonation  products 
=  volume  of  the  gaseous  detonation  products 

The  gas  pressure  varies  as  a  function  of  the  gas  volume.  If  we  assume 
a  polytropic  expansion  for  the  gas  then 

P  -  p  (Vf 

8  08  \\  j 

This  can  be  substituted  into  the  gas  expansion  work  term  to  give 


iP  dV  =  -  P  V 
I  g  g  og  og 


n  -n  . 

IV  c 

J  ^ 


The  work  done  on  the  canister  by  the  system  is  given  by 
W  =  ip  dV 

P  J  y  P 

where 

P^  =  pressure  required  to  exceed  the  metal  yield  stress 

V  =  volume  of  the  canister 
P 

Substitution  of  relations  (3),  (A),  and  (5)  into  Eq .  2  gives  the 
differential  form  of  the  equation  of  ’-otion  for  Regime  I  which  is 


P  dV  +  P  V 
y  P  og  og 


V  v 

n  f  ®  -n  f  a 

1  v  dV  ■  h,v 
)  8  g  1  *  • 


Since  the  mass  of  the  liquid  is  constant  and  if  we  assume  the  pressure 
required  to  expand  the  metal  to  be  a  constant  the  equation  of  motion 
can  be  integrated  directly  to  give 
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-P  (V  -V  )  ,  . 

y  p  op  1-n 


P  V  .  m. 

+  -g&-S8  (V  1-n  -  V  1_n)  =  -z-  v 


og  2  a 


or 


2  2 
/  =  — 
a  mrt 


p  v 

-?g  -°g  (V  1_n  -  V  1_n)  -  P  (V  -V  ) 
1-n  g  og  y  p  op 


(6) 


Since  we  are  interested  in  the  speed  of  the  leading  edge  of  the  cloud 
rather  than  in  the  average  speed  of  the  cloud  as  given  by  the  above  rela¬ 
tion,  we  assume  a  volume  relation  for  the  system  (see  Fig.  15),  such  that 


2  2 
ttL  <r  -  r  ) 
a 


cons  t . 


where 


L  =  canister  length 

r  «*  outside  radius  of  liquid  shell 

r  *  radius  at  which  v  occurs 
a  a 

Therefore,  half  of  the  liquid  volume  is  at  a  greater  radius  than  r  ,  and 
half  is  at  a  lesser  ratius. 


Differentiating  with  respect  to  time  we  have 


2,1  (r 


dr 


a 


0 


or 


v 


a 


since 


v 


a 


dr 
_ a 

dt 


v 


dr 

dt 


Note  that  the  average  cloud  speed  is  greater  than  the  speed  of  the  edge 
of  the  cloud.  The  liquid  particles  nearer  the  center  of  the  cloud  must 
be  moving  faster  than  those  near  the  edge  because  the  radius  of  the  cloud 
is  increasing  but  the  liquid  volume  is  assumed  to  remain  constant.  This 
is  illustrated  in  Fig.  16. 
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By  substitution,  the  relation  between  the  average  cloud  speed  and 
the  speed  of  the  edge  of  the  cloud  can  be  expressed  in  terms  of  volumes 
(see  Fig.  15) . 


Equation  6  can  further  be  simplified  by  noting  that 

V  1-n  «  v  ^-n 

8  og 

so  that  the  former  term  can  be  neglected  at  t  +  At. 

o 

Fur thermore 

V  =  V  -  V. 
g  p  £ 

so  that  the  equation  of  motion  for  Regime  I  becomes  finally 

V. 


,  2V 
v2  =  — £ 


mr  V 

£  P 


P  V 

-I-  P  (V  -  V  ) 
y  op  p 


n-1 


(7) 


Yield  Pressure 


In  order  to  determine  the  yield  pressure  of  the  metal  (Py)  we  assume 
that  the  relation  for  the  hoop  stress  in  a  thin  wall  cylinder  is  appli¬ 
cable  : 

P  -  S  - 
y  y  r 

where 

S  =  stress  at  which  metal  yields 

y 

=  pressure  at  which  metal  yields 
r  =  canister  radius 
t  =  canister  thickness 

If  we  assume  the  following  values  for  a  typical  three-pound  device 
Sy  =  35,000  psi  t  =  0.050  inch  (in  groove)  r  =  2  inches 
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the  yield  pressure  is  calculated  to  be  875  psi  or  1.26  x  103  psf. 
Calculation  of  Speed  in  Regime  I 

The  maximum  expansion  speed  in  Regime  I  can  be  calculated  by  the  use 
of  Eq.  7  since  we  know  or  can  calculate  all  the  required  values  except 
the  poly tropic  exponent,  n.  We  don’t  know  the  relationship  between  n  and 
the  gas  volume.  We  can,  however,  determine  an  average  value  of  n  by  noting 
that  the  ass'imed  value  of  n  was  2.0  in  the  Initial  Regime  and  that  its 
value  approaches  1.2  as  the  detonation  products  approach  perfect  gas  con¬ 
ditions.  We  will  assume  an  average  value  of  n  =  1.6. 

The  final  volume  can  be  obtained  by  assuming  a  linear  elongation  be- 
_fore  breakage  for  the  metal  of  the  canister.  The  final  volume^  is  given 
by 

2 

VT  —  V  £  £  =  maximum  elongation 

Ip  op 

Three-Pound  Device 

The  values  required  for  substitution  in  Eq.  7  are  obtained  from 
Table  1: 

V£  “  VP£  *  5-18  x  10_2  ft;3 

V  =  11.4  x  10"A  ft3 
°g 

V  =  V  +  V  =  5.29  x  10-2  ft3 
op  og  l 

P  =  4 . 6  x  10  2  psf 
og 

m^  =  0.093  slug 

P  =  1.26  x  105  psf 

y 

If  the  maximum  linear  elongation  is  assumed  to  be  1.2  then  the  final 
volume  is  given  by: 

VT  =  e2V  =  (1.2) 2  (5.29  x  10_2)  =  7.62  x  lo"2  ft2 

Ip  op 


3 

The  subscript  (I)  denotes  the  end  of  Regime  I. 
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These  values  can  be  substituted  in  Eq.  7  to  give 
2 


2(7.62  x  10~2)  -  5.18  x  1C~2 
-2. 


0.093(7.62  x  10  ) 


(4.6  x  IQ7) (11.4  x  10~4) 


1.6-1 


(1.26  x  103)  (5.29  x  10-2-  7.62  x  10”2) 
Simplification  then  gives 


v2  =  1.20  x  106  ft2/sec2 


or 


v  =  1100  ft /sec 

This  is  nearly  the  experimentally  determined  speed  of  the  fragments  as 
they  leave  the  canister  as  illustrated  by  Fig.  6.  The  value  of  v  is 
strongly  influenced  by  the  value  of  the  poly tropic  constant,  n.  Had  we 
chosen  values  of  1.5  or  1.7  (6%  change)  for  n,  the  calculated  speed  would 
have  varied  by  about  10%  which  would  still  be  a  reasonable  value  for  v. 

Eightv-Three  Pound  Device 

The  initial  values  required  for  insertion  in  the  equation  of  motion 
are  taken  from  Table  1  as  follows: 

Vi  =  m£/p£  =  1,43  ft;3 


°g 


V,  =  14  in3  =  8.11  x  10  3  ft3 

D 


V  =  V  +  V.  =  1,44  ft' 
op  og  l 


P  =  1.49  x  10  psf 

og 


m^  =2.59  slugs 

The  canister  yield  pressure  (Py)  is  calculated  from  the  equation 
for  hoop  stress  and  has  a  value  of  370  psi  or  5.33  x  1C ^  psf  based  on 
the  following  canister  characteristics: 


Yield  stress: 
Thickness  of  metal: 
Canister  radius: 


50,000  psi 
0.050  inches 
6.75  inches 


As  will  be  shown  later,  the  maximum  linear  elongation  prior  to  rup¬ 
ture  of  the  metal  in  the  canister  must  be  about  1.025  in  order  to  satisfy 
the  experimental  data  so  that 
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VT  =  c2v  »•-  (1.025) 2  (1.44)  =  1.513  ft3 
Ip  op 

Substitution  of  the  above  values  Into  Eq.  7  gives 

2  2(1.51)  -  1.43  [(1.49  x  108)(8.1  x  10~3) 

v  ”  2.59(1.51)  1  1.6-1 

5.33  x  10^(1 .44  -  1.51)] 

which,  on  simplification,  gives 

2  5  2  2 

v  =  8. 18  x  10  ft  /sec 


or 

v  *  904  ft/sec 

This  value  is  about  200  ft/sec  slower  than  that  calculated  for  the  thr.e- 
pound  device,  the  difference  being  primarily  a  result  of  the  lower  bur¬ 
ster  charge-to-liquid  ratio  used  for  the  83-pound  device  and  the  lower 
value  of  elongation  which  was  assumed  for  the  83-pound  device. 

Gas  Pressure  at  the  End  of  Regime  I 

As  a  further  check  we  can  calculate  the  gas  pressure  of  the  detona¬ 
tion  products  at  the  end  of  Regime  I  by  the  use  of  the  polytropic  equation  j 

for  the  gas  j 


Three-Pound  Device 

Substitution  of  the  appropriate  values  for  the  pressure  and  volume 

gives 


Ft  «  3.42  x  105  psf 

Ig 


or 

PT  -  2370  psi 

lg 


| 

A 

l 


j 


1 
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This  value  of  the  pressure  is  in  the  right  range  to  give  the  observed 
speed  of  the  squirt  spikes  as  shown  in  the  next  section. 

Eighty-Three  Pound  Device 

The  value  for  the  pressure  at  the  end  of  Regime  I  is  given  by  the 
polytropic  equation  where 


-  V,  =  1.513 


1.430  =  0.083  ft' 


PT  =  1.49  x  10 

Ig 


8.10  x  10 
0  .083 


-3  \  1.6 


PT  ■  3,6  x  10  psf 

Ig 

PT  =  25,000  psi 

Ig 

This  is  considerably  more  than  the  value  calculated  for  the  three-pound 
device.  If  we  refer  to  the  experimental  radius  versus  time  curves  for 
the  two  devices  (Fig.  8  and  12)  it  can  be  seen  that  the  initial  squirt 
jet  speed  for  the  83-pound  device  appears  to  be  much  greater  than  the 
three-pound  device.  In  fact,  the  initial  squirt  jet  speed  for  the  83- 
pound  device  appears  to  be  almost  2000  ft/sec  as  compared  to  some 
1225  ft/sec  for  the  smaller  device.  In  order  that  this  might  be  so,  the 
canister  pressure  at  the  end  of  Regime  I  for  the  83-pound  device  would 
have  to  be  on  the  order  of  25,000  psi  as  will  be  shown  later.  This 
could  be  so  if  the  canister  ruptured  sufficiently  early,  since  the 
internal  pressure  would  be  much  greater  under  those  conditions  at  the 
time  of  rupture. 

Squirt  Jet  Speed 

The  initial  speed  of  the  squirt  jets  can  be  approximated  from  the 
orifice  equation: 


v  =  C 
s  d 


where 


=  discharge  coefficient  (empirical) 
p.  =  liquid  density  (slugs  per  cubic  foot) 


A?  =  pressure  drop  across  the  orifice  (lb/ft  ) 
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For  this  particular  case  the  pressure  drop  can  be  assumed  equal  to  the 
pressure  of  the  detonation  products  minus  the  atmospheric  pressure  at 
the  time  of  rupture. 

Three-Pound  Canister 


The  stagnation  pressure  of  the  air  outside  the  canister  is  approxi¬ 
mately  two  atmospheres  or  2000  psf  at  the  time  of  rupture.  The  pressure 
drop  is  thus  given  by 

AP  «  PT  -  P  *=  3.42  x  105  -  0.02  x  105  =  3.40  x  105  psf 
Ig  a 

If  the  discharge  coefficient  is  assumed  to  be  0.5  (which  is  only  a  guess) 
and  the  liquid  density  is  1.80  slugs/ft-*  then 

I  - 

n  .  n  / 2(3.40  x  10s) 

vs  °'5  V L80 

v  ■  300  ft/sec 
s 

The  speed  of  the  jet  is  then  the  squirt  speed  plus  the  fragment  or  cloud 
speed 

v  =>  v  +  v  =  300  +  1100 
J  s 

Vj  •  1400  ft /sec 

This  is  to  be  compared  to  the  measured  initial  speed  of  the  squirt  3pikes 
which  is  about  1225  ft/sec  (see  Fig.  8) . 

Eight-Three  Pound  Device 

The  lack  of  available  information  for  this  device  makes  it  difficult 
to  analyse  the  squirt  jet  speed.  The  cloud  speed  and  the  pressure  at  the 
end  of  Regime  I  which  we  obtained  in  the  earlier  section  do  not  agree 
with  the  experimental  data  which  v;e  have  available.  But,  if  we  assume 
that  the  experimental  data  is  correct  and  assume  a  squirt  spike  speed  of 
about  2000  ft/sec,  the  pressure  drop  required  to  give  this  speed  can  be 
determined  from  the  orifice  equation.  If  we  assume  the  speed  of  the 
fragments  to  be  1000  ft/sec  the  speed  of  the  jet  emerging  from  the  can¬ 
ister  is 

2000  -  1000  »  1000  ft/sec 
The  orifice  equation  is  arranged  in  the  form 
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Substitution  of  the  appropriate  values  gives 


/ iooo\  2  lji 

=  \°‘5  /  2 


25,000  psi 


=  3,600,000  psf 


Elongation  of  the  Canister 

We  can  determine  the  maximum  canister  elongation  which  would  give 
the  above  pressure  at  the  time  of  rupture.  Knowing  the  initial  and  final 
pressures  and  the  initial  volume  we  can  calculate  the  final  volume  by  the 
polytropic  equation. 


VIg  " 


Substituting  the  appropriate  values  we  get 


1.49  x  10 
3.6  x  10£ 


8  \  1/1.6 


8.1  x  10-3  =  (41. 4)°’625  (8.1  x  10'3) 


VT  »  8.3  x  10"2  ft3 
Ig 

The  elongation  of  the  canister  is  given  by  the  equation 


where 


V  7  +  V. 

Ig  £ 


8.3  x  10-2  +  1.43  =  1.51  ft3 


=  1.44  ft“ 


so  that 


/1.05 


1.025 


This  suggests  that  eithet  the  canister  is  rupturing  almost  instantaneously 
without  stretching  appreciably  (possibly  as  a  result  of  impulsive  loading) 
or  that  the  squirt  speed  measured  from  Fig.  12  is  wrong.  The  character¬ 
istics  of  aluminum  and  steel  under  impulsive  loads  provide  evidence,  however, 
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that  the  steel  canister  might  well  rupture  earlier  than  the  aluminum 
canister.  Reference  5,  for  example  shews  that  the  critical  impact 
velocity  for  aluminum  is  about  twice  that  for  steel: 


Material 


24 ST  Aluminum 


1020  Steel 


Critical  Impact  Velocit 


200+  ft/aec 
100  ft/9ec 


The  critical  impact  velocity  is  the  speed  at  which  a  material  shatters 
at  the  point  of  impact.  As  such  it  is  a  measure  of  the  resistance  of  a 
metal  to  shattering  under  an  impulsive  load  such  as  a  detonation  wave. 

The  lower  value  of  critical  impact  velocity  for  steel  indicates  that  it 
is  more  likely  to  shatter  under  a  given  impulsive  load  than  is  aluminum. 

If  the  canister  of  the  83-pound  device  does  rupture  at  an  elongation 
of  only  1.025  then  the  fragment  speed  at  the  end  of  Regime  I  will  be 
lower  than  for  the  three-pound  device.  We  know  that  the  fragments  of 
the  83-pound  device  do  travel  considerably  further  than  the  fragments 
for  the  three-pound  device.  This  can  be  explained  by  the  greater  mass 
per  unit  area  of  the  steel  fragments  of  the  83-pound  device  as  compared 
with  the  aluminum  fragments  of  the  three-pound  device.  It  might  be  in¬ 
formative  to  determine  experimentally  the  distance  traveled  by  fragments 
from,  say,  a  steel-walled  three-pound  canister  or  an  aluminum-walled 
83-pound  canister. 


REGIME  II 

During  Regime  II  we  assume  that  a  "fluid  piston"  composed  of  the 
still  expanding  detonation  products,  the  liquid  shell  surrounding  the 
gaseous  products,  and  the  metal  canister  fragments  governs  the  outward 
expansion.  The  fluid  piston  is  assumed  to  be  a  coherent  mass  and  the 
liquid  surrounding  the  deto’  tion  products  is  assumed  not  to  have  broken 
into  discrete  droplets,  but  is  able  to  contain  the  detonation  gases 
within  it.  The  average  speed  of  the  expanding  cloud  is  increasing  dur¬ 
ing  this  regime  because  the  net  expansion  force  is  positive.  The  gas 
pressure  is  greater  than  the  atmospheric  pressure  and  there  is  an  out¬ 
ward  acceleration  of  the  liquid. 

During  Regime  II  the  cloud  is  divided  into  two  distinct  parts  which 
expand  at  different  rates.  They  are: 

1.  The  squirt  spikes  are  assumed  to  instantaneously  reach  their 
maximum  speed  at  the  beginning  of  Regime  II  and  then  to  decelerate  as 
a  result  of  aerodynamic  drag  in  the  disturbed  atmosphere  behind  the  pro¬ 
pagating  shock  wave.  The  liquid  in  the  squirt  spikes  is  assumed  initi¬ 
ally  to  be  in  the  form  of  a  jet  of  liquid  being  forced  outward  through 
an  orifice  formed  by  the  separating  walls  of  the  canister. 


36 


I  ir umiir *»  i»irwiijjrrr»T«ivi*tf«^r,inr^,T^«|i^ir'W^,,''l|^^^W?nTWW^!7^r«ll^ll»iPI!iuiifilf«|ii'11  !""wn*TjW*www,'wrr7 


NWC  TP  4702 


2.  The  main  Body  of  the  cloud  is  assumed  to  be  the  remainder  of  the 
liquid  between  the  fragments  and  the  squirt  spikes.  Initially  it  is 
assumed  to  be  a  coherent  mass  which  contains  the  detonation  products 
within  it.  Gradually  It  breaks  up  into  droplets,  but  in  this  regime 
the  droplets  are  packed  sufficiently  close  together  so  that  they  entrain 
any  air  between  them  and  continue  to  act  as  a  coherent  mass. 

In  this  regime  we  assume  the  speed  of  the  edge  of  the  main  body  to 
be  the  speed  the  main  body  would  achieve  if  it  were  not  restrained  by 
the  fragments .  The  main  body  of  the  cloud  is  forced  outward  by  the  pres¬ 
sure  of  the  detonation  products  and  in  turn  the  main  body  presses  outward 
against  the  fragments  and  the  atmosphere.  The  main  body  also  presses  out¬ 
ward  against  the  droplets  comprising  the  squirt  spikes,  an  effect  which  is 
assumed  to  be  negligible. 

We  assume  that  the  entire  mass  of  the  main  body  moves  uniformly  out¬ 
ward  and  that  any  atmospheric  air  which  is  entrained  moves  outward  at  the 
same  speed  as  the  liquid.  We  further  assume  that  the  volume  of  the  cloud 
is  the  same  as  the  original  liquid  volume^  and  that  the  cloud  expands  only 
in  the  radial  direction  (i.e.,  the  height  of  the  cloud  remains  constant  and 
equal  to  the  original  canister  length). 

At  the  beginning  of  Regime  II,  when  the  canister  ruptures,  the  strips 
move  outward  at  approximately  the  same  initial  speed  as  the  main  body  of 
the  cloud.  We  assume  that  during  Regime  II  the  fragments  accelerate  more 
slowly  than  the  droplets  comprising  the  main  body  of  the  cloud  since  they 
are  more  massive.  Thus  a  portion  of  the  cloud  reaches  a  greater  maximum 
speed  than  the  fragments.  A  problem  arises  here  since  Regime  II  extends 
outward  to  only  two  or  three  canister  radii.  Therefore,  a  considerable 
portion  of  the  main  body  may  be  trapped  behind  the  fragments  and  cannot 
accelerate  more  rapidly  than  the  fragments.  We  will  assume  that  the 
speed  of  the  main  body  is  the  speed  which  the  main  body  would  achieve  if 
it  were  not  restrained  by  the  fragments.  Note  that  even  though  the  speed 
of  the  leading  edge  of  the  cloud  might  be  restrained  from  increasing 
markedly  by  the  fragments,  the  average  speed  of  the  cloud  can  still  in¬ 
crease  markedly  because  the  droplets  closer  to  the  center  must  travel 
faster  than  the  droplets  nearer  the  outer  edge  in  order  to  satisfy  the 
condition  that  the  liquid  volume  remain  constant. 

Regime  II  begins  at  the  point  of  canister  rupture  and  continues  to 
the  point  at  which  the  main  body  of  the  cloud  reaches  its  average  maxi¬ 
mum  speed. 

In  this  regime  the  canister  has  ruptured  so  the  canister  expansion 
work  term  disappears  from  the  equation  of  motion.  Furthermore,  the 
liquid  which  has  squirted  out  through  the  breaks  in  the  canister  wall 


in  order  to  maintain  a  constant  volume  liquid  shell,  the  thickness 
of  the  liquid  shell  decreases  as  the  radius  of  the  system  increases. 
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must  be  accounted  for  In  terms  of  lost  mass  and  lost  energy  since  we 
assume  these  squirt  spikes  to  be  independent  of  the  main  body  of  the 
cloud.  A  control  volume  similar  to  that  postulated  for  Regime  I  is 
used. 

Equation  of  Motion 

The  energy  equation  for  a  non-flow  system,  which  is  again  used  as 
the  basis  for  the  equation  of  motion,  can  be  written 

fp  dV  -  Cp  dV  m.  (vdv  +  E 

Jgg  J  as  p  £pj  s 

where 

P  *»  atmospheric  pressure  at  the  face  of  the  piston^ 


*p 


energy  lost  to  the  squirt  spikes 
mass  of  liquid  in  fluid  piston 


We  integrate  from  the  point  of  canister  rupture  (j)  to  the  point  of 
maximum  speed  and  the  end  of  Regime  II  (jj)  .  After  substitution  of  the 
polytropic  relation 


P  =  PT 
g  Ig 


we  get  the  equation 
P 


L  _  v  _  p  (v  _v  )_!!&£  L 

1-n  yg  Ig  J  asV  p  Ip;  2  \a 


la 


+  E 


The  speed  va  is  the  average  speed  cf  the  liquid  mass  at  any  volume,  V. 
Solving  for  the  square  of  the  speed,  v  we  get 


P  V 

(V  1_n_  v  1_n)  _  p  (v  -V_  )  -  E 
1-n  g  Ig  as  p  Ip  s 


+  v 


la 


Pas  is  a  function  of  the  cloud 
pressure  at  the  face  of  the  piston, 
in  this  Regime  we  will  assume  P  to 
gration. 


speed  since  it  is  the  stagnation 
Since  the  speed  is  fairly  constant 
be  constant  to  simplify  the  inte- 
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Since  we  want  the  speed  c£  the  edge  of  the  cloud  rather  than  the  average 
speed  we  substitute 


v 


to  get  the  equation  of  motion  for  Regime  II: 

r  n 

i&_  (v  1"n-  vT  1-n) 

n  g  Ig 

+  vig/ 

Energy  Lost  to  Squirt  Spikes 

The  energy  lost  to  the  squirt  spikes  is  a  function  of  the  squirt 
spike  speed  and  the  mass  of  the  liquid  in  the  squirt  spikes. 

Three-Pound  Device 


+  V 

2 _ g_ 


2-p 


P  _  ' 


P  (V  -V-  )  -  E 
as  p  Ip  s 


The  speed  of  the  squirt  spikes,  v  j  ,  was  calculated  in  an  earlier 
Section  as  1400  f c/sec.  We  don't  know  the  mass  of  liquid  in  the  squirt 
spikes,  but  estimate  it  to  be  about  10%  of  the  total  liquid  mass,  or 

=  O.lm^  =  0.1(0.093  slug)  =  0.0093  slug 


The  energy  lost  to  the  squirt  spikes  is  then 


Is  ,  2  2, 

T  (Vj  -  v  > 


where 

Vj  =  squirt  spike  speed 
v  =  main  body  speed 
substituting  numerical  values  we  get 


E  =  ■°-°°9'3  (14002  -  11002) 

8  <- 
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E  -  3500  ft-lb 
s 

Eightv-Three  Pound  Device 

If  we  again  say  that  10%  of  the  total  liquid  mass  goes  into  the 
squirt  spikes  the  energy  lost  is  given  by 

E  =  (20002  -  10002) 

9  l 

E  =  3.9  x  105  ft-lb 

s 

This  is  an  awful  lot  of  energy  to  be  lost,  but  it  follows  naturally  from 
the  experimental  data  which,  if  correct,  predicts  a  squirt  spike  speed 
of  2000  ft/sec  for  the  83-pound  device. 

Volume  at  Which  Maximum  Expansion  Speed  Occurs 

The  maximum  expansion  speed  for  the  fluid  piston  occurs  when  the 
pressure  of  the  detonation  products  inside  the  piston  equals  the  stagna¬ 
tion  pressure  of  the  atmosphere  outside  the  piston.  At  this  time  the  sum 
of  the  forces  is  equal  to  zero  and  the  acceleration  of  the  piston  is  equal 
to  zero.  The  poly  tropic  relation  for  Regime  II  can  be  used  to  find  the 
cloud  volume  at  which  the  expansion  speed  is  a  maximum.  At  that  point 
Pg  =  Pgg  and  we  consider  this  the  end  of  Regime  II.  The  poly tropic 
equation  is 


In  this  regime  the  expanding  detonation  products  are  approaching 
perfect  gas  conditions  so  that  the  exponent,  n,  will  be  approximately 
1.3. 


The  atmospheric  pressure  can  be  calculated  from  the  gas  dynamic 
equation  for  compressible  flow: 


<9) 


where 


a  =  sonic  speed 
K  *  gas  constant 
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To  find  the  atmospheric  stagnation  pressure  we  assume  the  values 
P&  =  ambient  pressure  **  2120  psf 
K  *1.4  for  air 
v  =  1300  ft/sec  (calculated) 
a  =  1140  ft/sec 

Substitution  into  the  above  equation  gives 

1.4 


P  =  2120 
as 


1  +  0.4/2 


/l3Qo\ 2  1.4-1 

\1140/ 


P  =  4770  psf 
as  r 


The  initial  volume  of  the  gaseous  detonation  products  is  given  by 

VT  =  VT  -  V  -  8.95  x  10~2  -  5.18  x  10-2 
Ig  IP  £ 

VT  -  3.77  x  10”2  ft3 

Ig 

and  the  initial  pressure  of  the  detonation  products  is  1.70  x  103  psf. 

Substitution  in  the  polytropic  equation  results  in 

t  /i  in  lr,5  \  1/1 . 3 
VTT  =  (3.77  x  10~2)  I1'70  X-  10, 

IIg  y4. 77  x  103 

VIIg  =  °-589  ffc3 


The  volume  of  the  fluid  piston  or  main  body  of  the  cloud  is  then 

0.589  +  C .047 

x, 

VTT  =  0.636  ft3 
up 


VTT  =  VTT  +  0.9V. 
Up  IIg  l 


Note  that  we  have  accounted  for  the  loss  of  iiquid  to  the  squirt 
spikes  (about  10%)  in  making  this  calculation. 
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The  radius  of  the  fluid  piston  or  main  body  of  the  cloud  at  this  time  is 
then 

■i  A  Up  _  /Q-636 

rxip  V  1,L  '  y  TT  X  0 .833 

rTT  "  0 .49  ft 

rip 

Note  that  the  cloud  radius  at  which  maximum  speed  occurs  is  little  more 
than  one  canister  diameter.  This  accounts  for  the  fact  thac  any  increase 
of  cloud  speed  is  difficult  to  detect  in  the  experimental  data. 


Eighty-Three  Pound  Device 

The  volume  at  which  maximum  expansion  speed  occurs  can  be  calculated 
for  the  eighty-three  pound  device  in  a  manner  similar  to  that  for  the 
three-pound  device.  If  we  assume  that  the  maximum  cloud  speed  will  be 
approximately  the  same  as  for  the  three-pound  device  the  dynamic  pres¬ 
sure  outside  the  cloud  will  be  4770  psf.  The  gas  volume  at  maximum 
speed  will  then  be 


A  \ 1/n 

VT,  »  VT  I—2— )  -  8.3  x  10 

118  18  \PIIg/ 

VTT  «  13.7  ft3 

Ug 


(3.6  x  106 
4.77  x  103 


1/1.3 


The  piston  volume  will  be 
V T T  =  VTT  +  0.9  V. 

iip  iig  l 

VTT  =  15.0  ft3 
lip 


13.7  +  0.9  (1.43) 


and  the  main  body  radius  will  be 


Up 


,  /VIIp  _  / 15.0 

Y  ttL  y  n  1.6 


rTT  =  1.73  ft 
Up 


Maximum  Cloud  Speed  In  Regime  II 


The  maximum  cloud  speed  can  be  obtained  by  substituting  the  volume 
at  which  maximum  speed  occurs  (VjIp)  in  the  equation  of  motion  for 
Regime  II  (Eq .  8),  and  substituting  the  other  appropriate  values. 
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Three-Pound  Device 

The  following  value  are  given: 


V£ 

(5.18  x  10-2)  0.9  - 

VIg 

= 

3.77  x  10~2  ft3 

VIP 

- 

8.95  x  10'2  ft3 

m£ 

s 

0.093  (0.9)  =  0.084 

PIg 

s 

1.70  x  105  psf 

2 

VI 

= 

1.25  x  10^  ft2/sec2 

V__ 

_ 

0.589  ft3 

lig 

v._ 

= 

0.636  ft3 

lip 

p 

as 

a 

4770  psf 

n 

= 

1.3 

E 

s 

2500  ft- lb 

,-2  ,3 


Substituting  in  Eq.  8  we  get 

-2 


4.66  x  10 

2 _ 

0.636 


+  0.589 


2 

0.084 


(1.70  x  10  ) (3,77  x  10 
1  -  1.3 


-2,1.3 


(0.589-0-3  -  0.0377-0*3)  -  (4.77  x  103) (0.636  -  0.0895)  -  2500 

(8.95  x  10~2) (1.25  x  102) 

4.66  x  10-2  .  .  „  ,„-2 

- 2 -  +3.77  x  10 


2 


'II 


.92  x  10^  ft'Vsec^ 


v  =  1380  ft/sec 


A  second  iteration  would  improve  the  value  for  the  maximum  speed  if  this 
value  were  substituted  back  into  Eq .  9  but  the  improvement  would  not  be 
significant . 

This  is  the  maximum  speed  of  the  edge  of  the  cloud  for  the  three 
pound  device  in  Regime  II,  and  occurs  at  the  end  of  the  regime.  From 
this  point,  the  speed  of  the  cloud  decreases. 

It  is  interesting  to  compare  the  increase  of  the  average  cloud  speed 
from  Regime  I  to  Regime  II  with  the  increase  of  the  cloud  edge  r-peed  from 
Regime  I  to  Regime  II.  The  average  cloud  speed  at  the  end  of  Regime  I  is 


al 


al 


1100 


8.95  x  10 


5.18  x  10 

2 


_  o 


+  3.77  X  10 


-2 


and  the  average  cloud  speed  at  the  end  of  Regime  II  is 

=  1380 


all 


0 .642 


4.66  x  10 


-2 


+  0.589 


v  TT  *  1415  ft/sec 
all 


Note  that  the  edge  of  the  cloud  increased  some  280  ft/sec  in  speed  where¬ 
as  the  average  speed  increases  only  115  ft/sec  as  shown  by  the  table 
below: 


Average  Speed  Cloud  Edge  Speed 
Regime  I  1300  1100 

Regime  II  1415  1380 

Eighty-Three  Pound  Device 

The  maximum  cloud  speed  at  the  end  of  Regime  II  can  be  calculated 
in  a  manner  similar  to  that  for  the  three-pound  device  by  substituting 
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the  appropriate  values  as  follows: 


£ 

1.29  ft3 

Ig 

= 

8.3  x  10"2 

Ip 

= 

1.51  ft3 

\ 

= 

2.32  slugs 

Ig 

= 

3.6  x  106 

r  2 

I 

= 

8.18  x  105 

Ilg 

C 

13.7  ft3 

up 

= 

15.0  ft3 

as 

t= 

4770  psf 

i 

=s 

1.3 

= 

3.9  x  105 

Substitution  in  Eq.  8  resulcs  in 
,  +13-7  (  „ 


[3.6  x  10  ) (8.3  x  10 
1  -  1.3 


(13.7' 


-  0.083~0,3)  -  (4.47  X  103)(15.0  -  1.51)  -  3 


1.51(8.18  x  10' 


.9  x  105 


+  0.083 


B  1,91  x  10^  ft2/s ec2 


v  °  1380  ft/sec 


REGIME  III 


In  Regime  III  the  cloud  will  begin  to  decelerate  since  the  pressure 
within  the  piston  is  less  than  the  pressure  outside  the  piston.  The  gas 
within  the  piston  will  continue  to  expand  and  its  pressure  will  decrease 
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until  some  limiting  value  of  over-expansion  is  reached  as  a  result  of 
leakage  from  the  outside,  at  which  time  the  pressure  will  begin  to  re¬ 
turn  to  atmospheric.  The  limiting  value  of  the  pressure  within  the 
piston  is  the  pressure  as  expressed  by  the  polytropic  relation  for  any 
given  expansion  ratio.  This  latter  pressure  could  never  be  reached  in 
an  actual  situation  because  the  piston  is  becoming  less  coherent  and 
is  probably  breaking  up  into  discrete  liquid  particles.  There  will  be 
a  leakage  of  air  into  the  center  of  the  piston  from  the  top,  bottom, 
and  through  the  walls  of  the  fluid  piston.  This  leakage  will  have  two 
effects : 

1.  The  overexpansion  process  will  be  less  pronounced  since  the 
pressure  differential  across  the  piston  will  be  reduced,  and 

2.  The  liquid  particles  will  be  moving  at  a  relative  speed  with 
respect  to  the  air  leaking  into  the  center  of  the  piston,  which  will 
cause  aerodynamic  drag.  In  order  to  simplify  the  analysis  we  are  forced 
to  neglect  this  drag  in  Regime  III.^ 

At  the  same  time  the  stagnation  pressure  outside  the  piston  will  de¬ 
crease,  although  at  a  lesser  rate,  because  the  speed  of  the  expanding 
piston  is  decreasing.  The  minimum  value  of  the  stagnation  pressure  is 
one  atmosphere  for  zero  cloud  speed. 

During  this  regime  the  fragments  will  begin  to  pull  ahead  of  the 
main  body  of  the  cloud  and  to  pass  the  squirt  spikes.  The  fragments  are 
never  immersed  in  the  cloud.  They  are  always  at  the  cloud  edge  for  their 
particular  angular  position.  Some  of  the  droplets  of  the  main  body  will 
have  moved  outward  further  between  the  fragments  during  Reg:  II  and 

the  squirt  spikes  will  have  been  at  greater  radius  earlier. 

It  is  during  this  regime  that  the  fragment  spikes  begin  to  form. 

Since  the  fragments  have  relatively  large  surface  areas  they  form 
wakes  behind  them.  As  the  fragments  move  outward  with  increasing  speed 
relative  to  the  main  body  of  the  cloud  they  entrain  some  of  the  droplets 
in  their  wakes  and  carry  them  outward  at  speeds  greater  than  the  remainder 
of  the  cloud.®  Those  droplets  form  the  fragment  spikes.  The  fragments 
gradually  leave  the  droplets  behind,  however,  and  the  spikes  then  de¬ 
celerate  rapidly  in  the  ambient  atmosphere. 


7The  apparent  effect  of  this  simplifying  assumption  is  only  specula¬ 
tive,  but  appears  negligible. 

O 

The  mass  loss  to  the  fragment  spikes  is  not  accounted  for  in  the 
remaining  analysis  since  it  is  essentially  negligible  in  comparison  to 
the  total  mass  for  these  specific  designs. 
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Equation  of  Motion 

The  equation  of  motion  will  have  the  same  general  form  as  that  for 
Regime  II.  The  only  forces  acting  on  the  piston  will  be  the  pressure 
force  from  the  detonation  products  pressing  outward  and  the  stagnation 
pressure  force  of  the  atmosphere  outside  the  piston. 

Although  the  cloud  is  probably  breaking  up  into  discrete  droplets 
we  continue  to  treat  the  cloud  as  a  fluid  piston  and  ignore  the  effects 
of  aerodynamic  drag  on  the  individual  droplets  as  a  result  of  leakage 
through  the  piston.  It  is  impossible  to  determine  how  serious  an  error 
this  introduces,  but  the  complexity  of  attempting  to  deal  with  an  equa¬ 
tion  of  motion  treating  both  the  fluid  mass  and  the  individual  droplets 
in  the  same  equation  makes  it  necessary.  The  analysis  is  semi -empirical 
in  this  regime  because  we  can't  determine  the  influence  of  such  factors 
as  leakage  of  air  from  the  outside  atmosphere  into  the  inside  of  the 
pis  ton . 

Regime  III  is  probably  the  most  important  phase  of  the  expansion 
process  and  yet  is  the  one  in  which  we  are  least  able  to  predict  the 
results  analytically.  About  80%  of  the  expansion  of  the  cloud  occurs 
in  this  regime. 

In  the  case  of  the  three-pound  device  the  experimental  data  indicates 
that  the  relation  between  the  radius  of  the  cloud  and  the  time  seems  to 
change  at  a  time  of  about  5.5  msec  (see  Fig.  9).  From  Fig.  12  the  same 
phenomenon  seems  to  take  place  at  a  time  of  about  30  msec  in  the  case  of 
the  83-pound  device.  This  appears  to  be  the  point  where  the  fluid  piston 
effect  breaks  down  and  some  other  effect  becomes  predominant.  We  will 
assume  that  Regime  III  ends  at  this  point  and  that  another  regime  with 
a  different  equation  of  motion  begins.  The  apparent  dividing  line  be¬ 
tween  Regime  III  and  Regime  IV  occurs  at  cloud  volume  of  about  55  ft^ 
in  the  case  of  the  three-pound  device  and  a  volume  of  about  2600  ft^  in 
the  case  of  the  83-pound  device  (see  Fig.  11  and  13). 

The  equation  of  motion  for  Regime  III  is  obtained  from  the  energy 
equation 


r 

£ 
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fc 

i 
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i 

J 

1 
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P  dV  -  P  dV  ■  m„vdv 
g  g  as  p  l 


(10) 


Since  the  pressure  outside  the  piston  i3  the  stagnation  pressure  we  re¬ 
write 


P 

as 


P 

a 


(1  + 


K-i 


2 

v 

2 


a 


K/K-l 


The  presence  of  the  speed  term  in  the  equation  for  stagnation  pressure 
and  also  in  the  kinetic  energy  relationship  makes  the  equation  of  motion 
quite  difficult  to  integrate  in  that  form.  It  would  be  possible  to  re¬ 
sort  to  numerical  techniques  to  integrate  the  equation,  but  since,  as 
will  be  shown  later,  an  empirical  relation  is  necessary  to  describe  the 
gas  pressure  (P  )  versus  volume  relationship,  it  did  not  seem  practical 
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to  try  to  preserve  the  analytical  form  of  the  stagnation  pressure  equa¬ 
tion.  It  was  more  practical  to  resort  to  an  empirical  relation  for  that 
term  also  and  preserve  the  capability  of  integrating  the  equation  of  motion. 

Based  on  that  reasoning,  it  was  decided  to  derive  an  empirical  re¬ 
lationship  between  atmospheric  pressure  and  volume  using  the  experi¬ 
mental  data  for  the  three-pound  and  83-pound  devices  to  determine  the 
atmospheric  pressure-volume  relationships. 

Furthermore,  the  shape  of  the  experimental  cloud  speed  versus  volume 
curve  (as  illustrated  by  Fig.  13)  suggests  that  the  rate  of  decrease  of 
the  speed  with  respect  to  the  volume  must  decrease  with  an  increase  in 
volume.  Hence,  the  pressure  pushing  inward  on  the  piston  must  become 
less  as  expansion  proceeds.  Because  the  pressure  of  the  detonation  pro¬ 
ducts  pushing  outward  is  by  definition  the  same  as  the  stagnation  pres¬ 
sure  of  the  atmosphere  at  the  beginning  of  this  Regime,  the  net  force 
must  increase  rather  than  decrease  (at  least  during  the  first  portion  of 
the  Regime) .  It  was  further  assumed  that  the  pressure-volume  relationships 
are  as  shown  schematically  in  Fig.  17.  During  the  first  portion  of  Regime 
III  (identified  as  Regime  IIIA)  the  net  retarding  force  increases  because 
the  gas  pressure  inside  the  piston  decreases  more  rapidly  than  the  stagna¬ 
tion  pressure  outside  the  piston  decreases.  The  gas  pressure  inside  the 
piston  drops  to  atmospheric  pressure  and  then  overexpands  slightly  because 
of  inertia  effects  before  starting  to  rise  back  to  atmospheric  pressure. 

The  atmospheric  stagnation  pressure,  on  the  other  hand,  decreases  smoothly 
to  some  pressure  greater  than  ambient  atmospheric  pressure  at  the  end  of 
the  Regime.  The  second  part  of  Regime  III  we  call  Regime  II IB.  The 
empirical  equation  of  motion  will  be  different  for  each  of  these  sub¬ 
regimes  and  each  will  be  treated  separately. 

We  can  only  assume  a  value  for  the  terminal  gas  pressure  inside  the 
piston  at  the  end  of  Regime  IIIB  since  we  have  no  experimental  information. 
A  reasonable  assumption  is  one  atmosphere  since  leakage  of  air  through 
the  top  and  bottom  of  the  piston  would  tend  to  equalize  the  inside  and 
outside  pressures.  The  stagnation  pressure  outside  the  piston  does  have 
a  value  greater  than  one  atmosphere,  however,  since  the  piston  is  moving 
at  an  appreciable  speed  and  we  are  assuming  negligible  leakage  through 
the  face  of  Che  fluid  piston.  The  amount  of  overexpansion  must  also  be 
assumed,  mainly  on  the  basis  of  forcing  the  cloud  speed-volume  relation¬ 
ship  to  have  the  desired  value. 

Three-Pound  Device 


In  order  to  determine  an  equation  of  motion  for  the  three-pound 
device  it  was  first  necessary  to  find'  an  empirical  relation  between 
the  stagnation  pressure  of  the  atmosphere  and  the  volume  of  the  cloud. 
This  was  done  by  using  the  experimental  data  of  Fig.  11.  By  the  use  of 
the  cloud  speeds  obtained  from  Fig.  11  it  was  possible  to  calculate  the 
stagnation  pressure  as  a  function  of  volume  as  shown  in  Fig.  18.  The 
resulting  curve  was  then  fitted  by  an  empirical  equation  of  the  form 


48 


5000 


FIG.  18.  Dynamic  Pressure  as  a  Function  of  Cloud  Volume  Three-Pound  Device 
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A  +  BV 


A,  B  &  C  are  constant 


and  the  resulting  equation,  which  is  a  reasonably  good  fit,  is 

1/2 

P  =  4770  -  326  V  '  (11) 

as  p 

This  equation  is  valid  for  the  entire  Regime  III  for  the  three-pound  de¬ 
vice. 

Regime  IIIA 

The  relation  for  the  gas  pressure  inside  the  piston  has  different 
forms  for  Regimes  IIIA  and  IIIB.  For  Regime  IIIA  we  use  the  polytropic 
relation  which  is 


where  n  is  the  polytropic  exponent  if  there  are  no  leakage  effects.  In 
this  case  the  exponent  may  be  considered  a  modified  polytropic  constant 
which  also  accounts  for  the  effects  of  atmospheric  leakage  into  the  in¬ 
terior  of  the  piston.  The  effect  of  leakage  into  the  interior  of  the 
piston  is  similar  to  decreasing  the  value  of  n  since  the  leakage  would 
tend  to  lessen  the  drop  in  pressure  during  the  expansion  process.  This 
was  not  necessary  for  our  gross  calculations,  so  n  =  1.3  was  used. 

Substitution  of  Eq.  11  and  12  into  Eq.  10  results  in 


dV  -  \  (4770-326V  1/2)  dV  *  \  m„vdv 

g  J  P  P  J  * 


which  can  be  integrated  to  give  the  equation  of  motion  for  Regime  IIIA. 


P  V 

lie  Us 


(V  1_n-  VTT  1_n)  -  4770(V  -VTT  )  +  217 (V  3/2-VTT  ' 
g  Hg  p  Up  P  Up 


Substitution  of  the  appropriate  numbers  gives  the  numerical  form 
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v  -  23.8(-7970V  -  4770  V  +  217  V  +  12050)  +  1.90  x  10  (14) 
g  P  P 

To  find  the  piston  speed  at  the  end  of  Regime  II1A  we  select  an  assumed 
value  for  Pg  and  use  the  polytropic  relation  to  solve  for  Vg  at  the  end 
of  Regime  IIIA.  If  we  assume  that  the  gas  overexpands  to  1750  psf 


VTTT.  =  1.27  ft" 
IIIAg 

V  “V  +  V 

IIIAp  IIIAg  i 


1.32  ft' 


These  values  can  be  substituted  into  Eq.  14  to  give 


vTTT  2  =  1.88  x  10^  ft2/sec2 
IIIA 


IIIA 
Regime  IIIB 


1370  ft/sec 


In  this  sub-regime  the  relation  for  the  atmospheric  stagnation  pres 
sure  remains  the  same,  but  the  relation  for  the  gas  pressure  is  different 
We  assume  the  gas  pressure  increases  from  1750  psf  at  a  volume  of  1.32  ft 
to  2120  psf  at  55  ft^  according  to  the  relation 


3 


g 


1660  +  6.90V 


(15) 


Substitution  of  Eq.  15  into  the  equation  of  motion  gives' 


V 

S 


(1660  +6.90  V  )  dV  - 
P  P 


i  (4770  -  326  V  1/2)  dV 
P  P 

;HIAp 


IIAp 

fV 

1  m.vdv 

K 

IIIA 


which  can  be  integrated  to  give  the  equation  of  motion  for  Regime  IIIB. 
Making  the  appropiate  substitutions 


9 

We  assume  that  the  piston  volume  equals  the  gas  volume  in  this  sub¬ 
regime  since  their  difference  becomes  negligible  as  the  cloud  expands. 
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v2  *»  23.8  +  3.45  V  2  +  217  V  3/2  -  3110  V  +  370ol  +  1.88  x  106  (16) 

l  p  p  p  J 

To  find  the  piston  speed  at  the  end  of  Regime  IIIB  we  substitute  the  piston 
volume  at  the  end  of  Regime  III  which  is  assumed  to  be  55  ft^  and  solve 
for  the  speed. 

2  6  2  2 
v  »  0.26  x  10  ft  /sec 

v  =>  510  ft/sec 

3 

Using  the  same  methods  we  can  calculate  the  speed  for  volumes  from  1.32  ft 
to  55  ft3  and  plot  a  curve  of  speed  versus  volume.  This  curve  is  shown 
in  Fig.  11.  It  can  be  seen  that  the  curve  compares  well  with  the  experi¬ 
mental  data. 

Eighty-Three-Pound  Device 

The  oquatlon  of  motion  for  the  83-pound  device  is  similar  to  that 
for  the  thiee-pound  device.  The  schematic  relations  between  pressure 
and  volume  are  shown  in  Fig.  17,  which  is  valid  for  the  83-pound  device 
as  well  as  for  the  three-pound  device.  In  order  to  develop  an  equation 
of  motion  we  take  the  curve  for  the  main  body  speed  as  a  function  of  time 
and  determine  an  empirical  relation  between  the  atmospheric  stagnation 
pressure  and  volume  as  we  did  for  the  three-pound  device.  Using  the  data 
from  Fig.  13  for  main  body  speed  the  stagnation  pressure  cf  the  atmosphere 
ahead  of  the  cloud  was  calculated  by  means  of  Eq.  9.  This  is  plotted  in 
Fig.  19.  An  empirical  equation  of  the  form 

P  =  13620  -  7140  V  °'06  (17) 

as  p 

was  then  fitted  to  the  experimental  data.  As  can  be  seen  from  Fig.  19, 
the  fit  is  quite  good.  Equation  17  is  valid  for  both  Regime  IIIA  and 
Regime  IIIB  for  the  83-pound  device. 

Regime  IIIA 

With  the  substitution  of  the  polytropic  equation  for  the  variation 
of  the  gas  pressure  in  Regime  IIIA  and  Eq.  17  the  equation  of  motion  be¬ 
comes 
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lie  Ilf 


(V  1_n-  V  1 
g  Hg 


)  -  13620  (V  -  VTT  ) 
P  up 


+  6730  (V^06  -  VIIpL°6)  f  vn2  (15 

To  find  the  cloud  speed  at  the  end  of  Regime  IIIA  we  need  the  values  of 

Vg  and  Vp  at  the  end  of  Regime  IIIA.  If  we  assume  a  value  of  n  *  1.3 

and  if  we  allow  the  gas  to  again  expand  to  1750  psf  as  for  the  three- 

pound  device  we  can  calculate  V  TT  by  the  polytropic  equation  to  get 

lIIAg  • 


VTTTA  =  31.7  ff 
IIIAg 


V  =  V  +  V„ 

IIIAp  IIIAg  £ 


31.7  +  1.3  =  33.0  ft' 


Substitution  of  the  above  values  into  Eq .  18  along  with  the  known  values 
results  in 


(31.7 


13.7  0>3) 


-  13620  (33.0  -  15.0)  +  6730  (33. 01'06-  15. 01-06)  +  1.91  x  106 

v  2  =  1.88  x  106  ft2/sec2 
IIIA 

v.. .  i  =  1370  ft/sec 
IIIA 

Rexi'ne  IIIB 

In  Regime  IIIB  we  can  use  the  stagnation  pressure  relation  given  by 
Eq.  17,  but  a  new  relation  for  the  gas  pressure  inside  the  cloud  is  re¬ 
quired.  If  we  assume  that  the  gas  pressure  increases  from  1750  psf  at 
the  end  of  Regime  IIIA  to  atmospheric  pressure  (2120  psf)  at  the  end  of 
Regime  IIIB  the  linear  gas  pressure  versus  volume  relation  becomes 

P  =  1740  +  0.144V  (19) 

This  can  be  substituted  along  with  Eq.  17  into  the  equation  of  motion  to 
give 
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^  (1740  -  0 . 144V) dV  - 
V 

IIIAp 

or,  upon  integration 


(13620  -  7140V0 •06)dV  =  \  m£vdv 


IIIAp 


IIIA 


11880  (V  -  VTTTA  )  +  0.072  (V  -  V 


1  HA 

+  6730  (V  -  VTTT. 

p  IIIAp 


p  IIIAp 
1.06 


IIIAp 


) 


+  v 


IIIA 


(20) 


_Io  find  the  speed  at  the  end  of  Regime  IIIB  we  substitute  the  appropriate 
values  for  volume  and  speed  to  get 


2 

v _ 

2 

^11880  (2600 

IIIB 

2.32 

1.06 

,1.06,1 

+  6730 

(2600  ,D  - 

33  )j  + 

2 

U 

,  2  ,  2 

VIIIB 

■=  7.0  x  10 

ft  /sec 

VIIIB 

=  260  ft/sec 

-  33)  -  0.072  (26002  -  332) 
1.88  x  106 


The  relationship  between  speed  and  volume  calculated  from  Eq.  20  is  shown 
in  Fig.  13.  It  can  be  seen  that  the  agreement  is  quite  good. 

It  should  be  noted  that  the  agreement  between  actual  and  calculated 
values  is  not  so  Important  as  the  fact  that  the  same  equation  of  motion 
and  the  same  set  of  assumptions  were  used  tor  both  the  three-pound  and 
the  83-pound  analyses  and  that  agreement  was  reached  in  both  cases. 

Hence,  it  seems  that  the  present  approach  is  reacting  to  the  scaling 
factors  correctly. 


This  empirical  approach  was  used  only  as  a  preliminary  attempt.  It 
should  be  possible  to  insert  the  gas  dynamic  equation  in  the  equation  of 
motion  and  integrate  the  resulting  equation  by  numerical  techniques  to 
eliminate  the  empirical  relation  between  atmospheric  dynamic  pressure 
and  cloud  volume.  Considerable  work  will  be  required  to  determine  the 
actual  relations  between  internal  gas  pressure  of  the  detonation  pro¬ 
ducts  and  volume,  nowevet.  At  present  wc  have  no  information  concerning 
the  amount  of  leakage  into  the  interior  of  the  cloud,  the  amount  of  over 
expansion,  and  so  forth.  Until  we  can  get  this  information  there  is  no 
way  to  eliminate  the  empirical  relation  between  Internal  gas  pressure 
and  volume.  The  actual  relations  should  be  determined. 
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REGIME  IV 

In  the  previous  regimes  we  neglected  the  effect  of  aerodynamic  drag 
on  the  fluid  droplets  in  the  main  body  of  the  cloud  and  assumed  that  the 
fluid  piston  was  the  primary  factor  influencing  the  dispersion  process. 
In  Regime  IV  we  assume  that  the  fluid  piston  has  disintegrated  and  that- 
aerodynamic  drag  on  the  individual  droplets  completely  controls  the  dis¬ 
persion.  Neither  of  these  assumptions  is  correct,  of  course.  Actually 
both  aerodynamic  drag  and  the  fluid  piston  affect  the  disperson  process 
simultaneously.  In  order  to  make  the  model  susceptible  to  analysis, 
however,  it  was  necessary  to  simplify  it  in  this  manner.  The  regime  is 
assumed  to  begin  when  the  total  volume  of  the  cloud  with  respect  to  the 
liquid  droplets  becomes  sufficiently  large  so  that  the  effects  of  atmo¬ 
spheric  entrainment  become  negligible  and  ends  when  cloud  motion  for  all 
practical  purposes  ceases.  In  the  two  sizes  of  devices  studied  there  is 
a  definite  breakover  point  in  the  radius-time  curve  which  seems  to  in¬ 
dicate  the  end  of  Regime  III  and  the  beginning  of  Regime  IV. 

In  order  to  simplify  the  problem  we  further  assume  that  the  pres¬ 
sure  of  the  expanding  detonation  products  is  equal  to  the  pressure  of 
the  atmosphere,  that  the  atmosphere  is  motionless,  and  that  the  only 
force  acting  on  the  droplet  is  aerodynamic  drag. 

Equation  of  Motion 

The  equation  of  motion  for  the  droplet  is 


IF  =  F,  =  m.a  =  -  C,A  v  2 

d  d  /  add 


where 


m^  =  mass  of  liquid  droplet 
a  =  droplet  deceleration  =  v  ■— 

v  =  droplet  speed 

R  =  droplet  distance  from  center  of  cloud 


=  atmospheric,  density 
C =  droplet  drag  coefficient 


droplet  frontal  area 


■»  droplet  speed  with  respect  to  the  atmosphere 
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Since  the  Reynolds  number  is  approximately  1000  the  above  equation  should 
be  a  good  approximation  to  the  actual  drag  relation. 

We  can  rewrite  the  equation  in  the  form 


PaCdAd 


dR 


since  the  atmospheric  speed  is  assumed  to  be  zero.  If  the  atmospheric 

speed  were  not  zero  v  would  not  be  the  same  as  v,. 

d 

The  ratio  o l  the  droplet  area  to  the  droplet  mass  can  be  written 
A  r2 

- _d  _  TT  Cd  =  3  - 

md  4/3irr  ,^p,  ^Pdrd 

d  d 

=  droplet  radius 
=  droplet  density 


The  equation  can  be  rewritten,  therefore,  in  the  form 


0.375  PaCd 
pdrd 


dR 


Except  for  the  droplet  density,  each  of  the  quantities  on  the  right-hand 
side  of  the  equation  is  a  function  of  the  distance,  R.  Since  the  atmo¬ 
spheric  density,  the  drag  coeff icient(6) ,  and  the  droplet  radius  all 
decrease  with  time  or  Reynolds  number,  we  can  assume,  as  a  first  approxi¬ 
mation,  that  their  quotient  is  constant.  Therefore  we  can  write: 


0.375 


P  C 

a  d 


Pdrd 


const.  =  K 


Thus 


KdR 
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Upon  integration  we  get 


-  in  v  +  £n  v  =  K  (R 
e  d  e  Hid 


RIII) 


In  order  to  make  the  above  equation  compatible  with  the  equation  of 
motion  for  the  other  three  regimes  we  use 


and  the  equation  becomes 


£n  v, 
e  d 


£n  V,TTJ 

e  1 1  Id 


*1  (V1/2  -  V 

irL  '  III 


1/2. 


=  exp 


£n  v 
e  Hid 


£1  (v1/2  -  v 

ttL  v  III 


1/2. 


(21) 


This  is  the  equation  of  motion  for  a  single  cloud  droplet  in  Regime  IV. 
By  substituting  appropriate  values  for  the  variables  we  can  make  it 
applicable  to  either  the  three-pound  or  83-pound  device. 

Three-Pound  Device 

The  constant  K  is 


K  = 


CaCdAd 

2  m , 


0.375 


Vd 

pdrd 


We  can  substitute  representative  values  as  follows: 
=  0.0024  slug/ft2 

C,  =  0.80 


Pd  =  1.30  slug/ft 

r.  ■  0.0002  ft 

d 

The  droplet  radius  is  representative  of  values  which  have  been  obtained 
by  Stanford  Research  Institute  in  their  studies  (Ref.  2)  whereas  the 
drag  coefficient  is  representative  of  the  drag  coefficient  to  be  expected 
of  a  spherical  droplet  in  this  Reynolds  number  regime.  The  constant,  K, 
therefore  becomes 
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K  = 


0.375(0.0024) (0.80) 


2.00 


(1. 80) (0 .0002) 

Upon  substitution  of  the  following  values 
L  =  0.834  ft 


Hid 


III 


510  ft/sec 
55  ft3 


the  equation  of  motion  becomes 


v,  =  exp 
d 


6.2  - 


1.24  (V 


1/2 


'.42)] 


(22) 


By  the  use  of  Eq.  22  we  can  solve  for  the  droplet  speed  for  various  values 
of  V.  For  instance,  when  V  equals  80  ft^  the  speed  is  given  by 

v^  =  exp  J^6 . 2  -  1.24  (8.94  -  7 . 4 2 ) J  =  75  ft/sec 


Fig.  11  shows  the  calculated  speed  versus  cloud  volume  curve  for  Regime 
IV.  Agreement  with  the  experimental  data  is  good. 

Maximum  Cloud  Radius 


The  maximum  cloud  radius  can  be  obtained  by  allowing  the  droplet 
speed  to  decay  to  a  negligible  value  and  then  solving  for  the  volume. 
For  instance,  if  we  assume  that  the  maximum  radius  is  reached  when 
v  =  5  ft/sec  or 


v  =  0 .01  v 

II  III 

We  can  solve  for  the  maximum  volume,  V 


IV’ 


1/2 


IV 


III 


S,n  100 
e 

1.24 


v 


Jin 


1/2 


Illd 

/IVd 


IT 

7TL 


=  1/2 


IV 


124  ft' 


60 


NWC  TP  4702 


The  maximum  cloud  radius  is  then 


IV 


IV 


ttL 


124 


(0.83) 


RIV  =  6.9  ft 


This  value  agrees  fairly  well  with  the  available  experimental  data  as 
shown  in  Fig.  10. 

Eighty-Three  Pound  Device 

The  constant,  K,  becomes 

K  =  0.33 


when  the  following  representative  values  are  substituted: 

vTTT,  =  260  ft/sec 

1  lid 

3 

V  °  2600  ft 

IITd 

L  -  1.60  ft 

p  *  0.0024  slugs/ft3 

a 

p_,  =  1.80  slugs/fc3 

d 

C ,  =  0.40 

d 

r .  *=  0.0006  ft 

d 


The  larger  droplet  size  with  its  resultant  lower  drag  coefficient  can  be 
justified  on  the  basis  of  the  lower  burster  charge  to  liquid  ratio  of 
the  83-pound  device  as  compared  to  the  three-pound  device.  The  lower 
energy  per  unit  mass  of  liquid  resulting  from  the  lower  burster  charge 
to  liquid  ratio  would  tend  to  result  in  larger  droplets.  Experimental 
evidence  has  confirmed  that  the  larger  device  produces  larger  droplets. 

The  equation  of  motion,  therefore,  becomes 

Vj  «=  exp  [5.56  -  0.147  (V1/2  -  51)1  (23) 

u  1  1 

for  the  83-pound  device.  Figure  13  illustrates  the  cloud  speed  calcul¬ 
ated  by  the  use  of  Eq.  23  compared  to  the  cloud  speed  taken  from  experi¬ 
mental  data. 
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Maximum  Cloud  Radius 


If  we  allow  the  droplet  size  to  decay  to  the  same  value  as  for  the 
three-pound  device  (5  ft/sec)  then 


v  *  0 .0  2v 
dIV  dill 


and  the  equation  becomes 

„  1/2  _ 


in  50 
e 


IV 


0.147 


+  51 


Vlv  -  5930  ft' 


The  maximum  cloud  radius  is  then  given  by 
RIV  “ 


IV 


59  30 
1.6  TT 

33  ft 


This  value  is  somewhat  greater  than  the  experimental  value  given  by 
Fig.  12,  but  is  within  10%  at  a  time  of  0.10  seconds. 


SUMMARY 


The  equations  of  motion  which  have  been  developed  in  this  report 
match  the  experimental  data  which  is  available  for  two  sizes  of  explosive 
dispersion  devices.  Those  devices  represent  the  largest  and  smallest  de¬ 
vices  for  which  we  have  experimental,  information,  and  hence  the  equations 
seem  to  be  reacting  correctly  to  the  scaling  requirements. 

The  equations  are  based  on  the  concep.  of  a  fluid  piston  in  the 
first  three  regimes  of  motion.  This  fluid  piston  is  essentially  a  gas 
bubble  surrounded  by  a  liquid  shell.  As  the  gas  bubble  expands  it  forces 
the  liquid  outward  to  form  the  characteristic  explosive  dispersion  cloud. 

In  Regime  I,  in  which  the  canister  expands  before  rupturing,  the 
maximum  cloud  speed  (which  is  the.  same  as  the  expansion  speed  of  the 
metal  canister),  is  calculated  to  be  1100  ft/sec  for  the  three-pound 
canister.  This  is  almost ' exactly  the  same  as  the  initial  speed  of  the 
fragments  from  the  three-pound  canister  as  measured  experimentally.  No 
data  is  available  concerning  the  speed  of  the  fragments  from  the  83- 
pound  device,  but  the  calculated  value  of  904  ft/sec  is  reasonable  in 
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view  of  the  lower  burster  charge-to-liquid  ratio  of  the  larger  device 
and  the  fact  that  the  initial  elongation  of  the  steel  canister  for  the 
larger  device  appears  to  be  less  than  that  for  the  aluminum  canister 
for  the  smaller  device.  Since  the  canister  walls  are  accelerating  in 
this  regime  the  sooner  they  rupture  the  lower  their  maximum  speed  will 
be.  It  would  be  helpful  to  get  some  experimental  information  on  the 
true  time  of  rupture  of  the  steel  canister  of  the  83-pound  device  to 
verify  this  assumption.  It  might  also  be  helpful  to  compare  the  time 
of  rupture  for  both  aluminum  and  steel  canisters  for  devices  of  the 
same  size. 

As  a  result  of  the  earlier  rupture  of  the  steel  canister  the  pres¬ 
sure  inside  the  canister  at  the  end  of  Regime  I  will  be  greater  for  the 
83-pound  device  than  for  the  three-pound  device.  This  greater  pressure 
would  account  for  the  much  greater  initial  cloud  speed  which  was  measured 
for  the  83-pound  device.  Since  the  only  portion  of  the  cloud  which  can 
be  measured  very  early  in  the  expansion  is  the  squirt  jet,  the  experi¬ 
mentally  determined  initial  cloud  speed  of  about  2000  ft/sec  for  the 
83-pound  device  must  be  the  squirt  jet  speed.  This  is  much  greater  than 
the  1340  ft/sec  measured  for  the  three-pound  device  and  would  require  an 
internal  pressure  of  some  25,000  psi. 

In  Regime  II  the  liquid  cloud  is  still  accelerating  because  the 
pressure  inside  the  canister  is  greater  than  the  pressure  outside  the 
canister.  The  maximum  cloud  speed  is  reached  when  the  pressures  inside 
and  outside  are  equalized.  That  is  the  end  of  Regime  II.  The  equation 
of  motion  predicts  that  the  maximum  cloud  speed  for  the  three-pound  de¬ 
vice  will  be  1380  ft/sec  and  that  it  occurs  at  a  cloud  radius  of  about 
0.5  ft.  The  maximum  speed  for  the  83-pound  device  is  also  about  1380  ft/ 
sec  and  occurs  at  a  cloud  radius  of  about  1.8  ft.  These  speeds  and  dis¬ 
tances  are  compatible  with  the  experimental  information.  Note  that  the 
radius  at  which  the  maximum  speed  occurs  is  quite  small  as  calculated 
from  the  equation  of  motion.  Analysis  of  the  experimental  information 
reveals  that  it  is  very  difficult  to  detect  an  increase  in  speed  of  the 
cloud  initially  because  it  occurs  so  close  to  the  point  of  detonation. 
Only  on  very  high-speed  cameras,  such  as  the  Cordin  camera,  can  the 
initial  acceleration  of  the  cloud  be  detected. 

In  Regime  III  the  cloud  begins  tc  decelerate.  This  regime  is  the 
one  in  which  most  of  the  cloud  expansion  occurs,  but  it  is  the  one  which 
we  are  least  able  to  handle  analytically  because  so  little  information 
on  the  actual  processes  within  the  cloud  is  available.  By  assuming 
semi-empirical  relations  based  upon  what  knowledge  we  'nave  of  the  physics 
of  the  cloud  expansion  we  were,  however,  able  to  develop  an  equation  of 
motion  which  gave  reasonable  results  for  both  the  small  three-pound 
canister  and  the  large  83-pound  canister. 
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We  need  considerably  more  information  on  the  expansion  dynamics, 
including  overexpansion,  leakage,  and  the  individual  droplet  effects, 
before  we  can  replace  the  empirical  relationships  with  theoretical 
relationships.  We  think,  however,  that  once  this  information  is  avail¬ 
able,  the  basic  equation  of  motion  will  be  satisfactory.  The  speed- 
distance  relationships  obtained  from  the  equation  cf  motion  match  those 
obtained  from  experiment  quite  closely.  It  should  be  realized,  however, 
that  even  the  experimental  speed-distance  relationships  are  only  approxi¬ 
mations.  It  is  quite  difficult  to  get  reliable  speed-distance  values  by 
graphical  differentiation  of  the  available  radius-time  information. 

The  speed  at  the  end  of  Regime  III  was  calculated  to  be  510  ft/sec 
for  the  three-pound  device  and  260  ft/sec  for  the  83-pound  device.  While 
these  are  values  quite  close  to  those  which  would  be  expected,  the  nature 
of  the  equation  is  such  that  very  small  variations  in  the  numerical 
values  assumed  for  the  equation  result  in  large  variations  in  the  final 
calculated  cloud  speeds.  We  aren't,  for  example,  justified  in  assuming 
that  the  actual  droplet  speed  for  the  three-pound  device  is  greater  than 
that  for  the  83-pound  device  at  the  end  of  this  regime  based  upon  the 
calculations.  The  uncertainty  is  such  that  the  speeds  could  well  be 
reversed.  Hence,  the  final  speeds  indicate  only  that  the  equations  of 
motion  are  able  to  produce  reasonable  relationships  between  speed  and 
distance  when  reasonable  numerical  values  are  used  in  the  equations. 

In  Regime  IV  we  shift  our  analysis  to  look  at  the  individual  drop¬ 
lets.  We  abandon  the  fluid  piston  concept  and  concentrate  on  aero¬ 
dynamic  drag  of  the  droplets  as  the  controlling  factor.  In  the  analysis 
we  assume  that  the  fluid  piston  disintegrates  instantaneously  at  the  end 
of  Regime  III  and  aerodynamic  drag  of  the  droplets  takes  over  immediately 
as  the  controlling  factor.  This  is  not  really  the  case.  Aerodynamic 
drag  on  the  droplets  is  probably  influencing  the  expansion  beforl  the 
fluid  piston  disintegrates  completely.  The  assumption  does  seem  to  re¬ 
flect  the  existence  of  a  breakover  point  in  the  experimental  radius¬ 
time  curves  at  which  the  drag  seems  to  increase  markedly.  The  droplets 
at  that  point  seem  to  slow  down  almost  instantaneously,  as  if  the  aero¬ 
dynamic  drag  were  suddenly  becoming  predominant.  This  sudden  decrease 
is  shown  in  the  experimental  data  of  Fig.  11  and  13  and  is  reflected  in 
the  calculated  values  for  the  equations  of  motion.  This  seems  to  in¬ 
dicate  that  the  actual  transition  from  the  fluid  piston  effect  to  the 
aerodynamic  drag  effect  might  be  quite  abrupt  in  reality  as  well  as  in 
theory. 

The  measured  diameter  of  the  cloud  at  maximum  expansion  is  less 
than  that  predicted  by  the  equation  of  motion  in  the  case  of  the  83- 
pound  device  and  greater  i'ldii  that  predicted  by  the  equation  for  the 
three-pound  device.  A  major  reason  for  the  difference  is  the  use  of 
a  droplet  speed  of  260  ft/sec  at  the  beginning  of  the  regime  tor  the 
83-pound  device  and  a  speed  of  510  ft/sec  for  the  three-pound  device  as 
calculated  from  Regime  III,  and  the  interdependence  of  the  drag  coeffi¬ 
cient  and  the  droplet  radius. 
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What  is  important,  is  that  in  both  cases  the  values  obtained  are 
reasonably  close  to  the  experimental  values.  Whenever  assumptions  were 
made  they  were  based  on  physical  evidence  or  reasonable  extrapolations 
from  physical  evidence.  It  was  not  necessary  to  introduce  any  unfounded 
"fudge  factors"  in  order  to  make  the  equations  of  motion  balance  for  both 
the  three-pound  and  83-pound  devices,  and  multiple  iterations  would  im¬ 
prove  the  agreement  with  experimentally  observed  motion. 


CONCLUSIONS 


The  major  conclusion  which  we  have  drawn  from  this  study  is  the 
importance  of  the  fluid  piston  in  determining  the  size  to  which  the  ex¬ 
plosively  dispersed  cloud  will  expand.  It  appears  that  the  cloud  expands 
to  some  critical  size  as  a  result  of  the  piston  effect  at  which  point 
the  piston  effect  deteriorates  rapidly  and  the  expansion  rate  of  the 
cloud  rapidly  approaches  zero. 

This  critical  size  is  directly  related  to  the  size  of  the  device, 
and  we  believe  that  it  is  a  direct  function  of  the  thickness  of  the  liquid 
column  surrounding  the  gaseous  detonation  products.  When  this  liquid 
shell  becomes  too  thin  to  effectively  prevent  the  leakage  of  air  in 
large  quantities  through  the  piston  into  its  interior,  the  drag  of  the 
atmospheric  air  rushing  past  the  droplets  from  the  outside  of  the  piston 
to  the  inside  because  of  the  pressure  gradient  between  the  inside  and 
outside  becomes  the  dominating  factor  and  the  droplets  slow  down  almost 
immediately . 

If  this  is  the  case,  then  the  factor  dominating  the  size  of  the 
expanding  cloud  will  be  the  initial  thickness  of  the  liquid  column 
surrounding  the  burster  charge,  which  is  effectively  the  diameter  of 
the  explosive  dispersion  device.  The  greater  the  diameter  of  the 
device  the  greater  will  be  the  ultimate  size  of  the  cloud. 

A  certain  amount  of  burster  explosive  will  be  required  to  disperse 
the  liquid,  but  the  device  does  not  seen;  to  be  highly  sensitive  to  the 
burster  charge-to-liquid  ratio.  The  three-pound  device  had  twice  the 
burster  explosive  for  every  pound  of  liquid  compared  to  the  83-pound 
device,  but  the  cloud  diameter  of  the  smaller  device  was  obviously  not 
twice  the  size  of  that  for  the  larger  device.  The  most  obvious  effects 
of  the  higher  burster  charge-to-liquid  ratio  will  be  to  slightly  increase 
the  speed  of  the  expansion  and  to  break  the  liquid  into  smaller  droplets. 
Smaller  droplets  might  delay  the  decay  of  the  piston  by  providing  a  less 
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in  the  aerodynamic  regime  of  flight. 
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Another  effect  of  the  larger  burster  charge  ratio  will  be  to  force 
more  liquid  into  the  squirt  spikes.  Since  the  material  in  the  squirt 
spikes  might  be  wasted  material,  any  factor  which  takes  liquid  out  of 
the  main  body  and  forces  it  into  the  squirt  spikes  might  be  undesirable. 

The  influence  of  the  canister  material  also  appears  to  be  secondary. 
The  steel  canister  appears  to  rupture  earlier  than  the  aluminum  canister, 
forcing  more  material  into  the  squirt  spikes .  The  energy  required  to 
burst  the  canister  in  either  case,  however,  is  much  less  than  that  avail¬ 
able  to  accelerate  the  liquid  in  the  main  body  so  that  the  effect  of  the 
canister  material  on  the  dispersion  of  the  main  body  of  the  cloud  is 
slight  in  terms  of  expansion  distance.  Effort  should  be  devoted  toward 
developing  a  material  which  will  minimize  the  amount  of  liquid  in  both 
the  squirt  spikes  and  the  fragment  spikes.  A  canister  material  which 
would  rupture  into  minute  fragments,  thus  eliminating  the  squirt  spikes 
and  the  fragment  spikes,  and  which  would  also  provide  additional  material 
to  reinforce  the  shell  of  the  fluid  piston,  might  well  improve  the  per¬ 
formance  of  the  device. 

Since  the  overexpansion  of  the  gases  within  the  interior  of  the 
piston  strongly  influences  the  expansion  of  the  piston  by  creating  a 
negative  pressure  gradient  it  might  be  possible  to  increase  the  expansion 
distance  by  allowing  leakage  into  the  interior  of  the  piston  from  the  top 
and  bottom  during  the  latter  stages  of  the  expansion.  This  would  elimi¬ 
nate  the  adverse  pressure  gradient  so  that  the  negative  expansion  work 
would  be  less,  and  would  also  decrease  the  tendency  of  the  atmospheric 
air  to  force  its  way  into  the  piston  through  the  liquid  shell,  thus 
prolonging  Regime  III  of  the  expansion  process. 

In  order  to  increase  the  size  of  the  expanding  cloud  in  the  aero¬ 
dynamic  regime  of  flight  two  possibilities  present  themselves.  The 
first  of  these  is  the  use  of  a  more  dense  liquid  and  the  second  is  the 
use  of  a  more  viscous  material.  Both  of  these  would  increase  the  travel 
of  the  droplets  in  the  aerodynamic  regime,  although  the  increase  would 
not  be  dramatic. 


66 


REFERENCES 


S 

1.  Stanford  Research  Institute.  Explosive  Dispersion  of  Volatile  s§ 

Liquids,  by  R.  C.  Robbins.  Menlo  Park,  California,  Project  PAU-  g 

4047  Technical  Report  No.  8,  22  January  1964.  CONFIDENTIAL.  M 

2.  U.  S.  Naval  Ordnance  Test  Station,  Inyokern.  Fuel-Air  Explosive  3 

Clouds,  Part  2,  Explosive  Dispersion  of  Volatile  Liquids,  by  R.  C.  f. 

Robbins.  China  Lake,  Calif.,  NOTS ,  October  1965.  (NavWeps  Report  jj 

8660,  Part  2),  CONFIDENTIAL.  1 


3.  U.  S.  Naval  Ordnance  Test  Station,  Inyokern.  Exploratory  Develop¬ 
ment  Problems  in  FAX  Weaponization ,  Transactions  of  the  Fourth 
Symposium  on  Warhead  Research,  by  Richard  J.  Zabelka,  Douglas  G. 
Kinney,  and  Lloyd  H.  Smith.  China  Lake,  Calif.,  NOTS,  1965. 

(TP  3984),  CONFIDENTIAL. 

4.  Cook,  M.  A.,  The  Science  of  High  Explosives,  Reinhold,  New  York, 

1958,  p.  284. 


5.  Reinhard,  John  S.  and  John  Pearson,  Behavior  of  Metals  Under  Impul¬ 
sive  Loads,  Dover,  New  York,  1965,  pp.  16-17. 

6.  Schlichting,  Hermann,  Boundary  Layer  Theory,  Fourth  Ed.,  McGraw-Hill, 
New  York,  1960,  p.  16. 


67 


UNCLASSIFIED 


$oc»iniy  Classification  _ _ 

DOCUMENT  CONTROL  DATA  ■  R  4  D 


;5rrunl)'  i)ji  nlita/un  of  title,  body  i >f  nbrtrari  and  inde  iiijj  annotation  n- u-f  by  m tered  when  the  overall  report  It  flottilltd) 


i  cR'Gi-na’inG  activity  Ct'oiporslc  4u»>iof; 

2a.  REPOST  SECURITY  CLASSIFICATION 

Naval  Weapons  Center 

UNCLASSIFIED 

China  Lake,  California  93555 

26.  C  R  O'JP 

I  i  HtPOAt  TITLE 

EXPLOSIVELY  DISPERSED  LIQUIDS.  PART  1,  DISPERSION  MODEL 


4  DESCRIPTIVE  notes  (Type  ot  teport  and  inf  fun  vr  dote  t) 


%  A(J  ThORiS!  (Firrt  name,  middle  initial,  liii  n*m#> 


Richard  J.  Zabelka  and  Lloyd  H.  Smith 


«  RtPOST  DATE  .  .  . 

December  1969 _ 

|4-  CON  TRACT  OR  GRANT  NO 

6.  PflOJEC  T  NO 

AirTask  A35-532-050/216-1/W11-62X 


d 

IU  OllTRlBuTlON  STATEMENT 


7 O.  TOTAL  NO  OF  PAGES  76.  NO  Of  REFS 

68  6 
a«.  ORIGINATOR'S  REPORT  NUUPERlSl 

NWC  TP  4702,  Part  1 


j  Ofc.  CTmIA  REPORT  NOiSl  ( Any  other  number?  mot  rr  .y  be  at  it grted 
th it  report) 


THU  OOCUMENT  U  SUSJECT  TO  S FECIAL  EXPORT  CONTROL!  AND  EACH  TRANSMITTAL  TO  FOREIGN 
GOVERNMENTS  OR  FOREIGN  NATIONALS  MAY  SE  MADE  ONLY  WITH  PRIOR  APPROVAL  OF  THE  NAVAL 
WEAPONS  CENTER. 


i.  iii  — i  tiiiM.,'  NO'n 


1-  iPOKSOn.lO  Mill,*,.  ACTIVITY 

Naval  Air  Systems  Command 
Washington,  D.  C.  20360 


A  preliminary  analysis  of  the  dispersion  process  of  liquid  droplets  into  an 
aerosol  cloud  is  presented.  This  assessment  of  the  complex  phenomena  is  not  all 
encompassing,  but  is  believed  to  shed  some  light  upon  the  behavior  of  the  various 
mechanisms  which  govern  the  expansion  of  explosively  dispersed  liquids.  This 
analysis  does  provide  a  basic  understanding,  and  provides  data  for  preliminary 
design  estimates. 


The  equations  presented  are  based  on  the  concept  of  a  fluid  piston  in  the  first 
three  of  four  regimes  of  motion-  This  fluid  piston  i6  essentially  a  gas  bubble 
surrounded  by  a  liquid  shell.  As  the  gas  bubble  expands  it  forces  the  liquid  out¬ 
ward  to  form  the  characteristic  explosive  dispersion  cloud.  The  fluid  piston  de¬ 
generates  into  discrete  droplets,  whose  motion  is  governed  by  aerodynamic  forces  in 
Regime  XV. 

These  equations  of  motion  were  developed  from  experimental  data  collected  with 
two  different  size  devices.  The  agreement  with  this  data  indicates  that  scaling 
requirements  are  satisfied. 


DD  ,?r.,1473  1PAGr- 

S/N  010  I • 807- 6fiG I 


_UMCLASSmED 

Security  Classification 


l'«r. s  ■li'iSikiw'  u- . . . .  i immmir  "i*  *  ♦  •wam . . -MittM****- 


Q 

OC 

< 

o 


o 

cc 

»- 

in 

CD 

< 


s 

3 


a.  B 
a>  x 
o  oc 


c 

o 

o 

<d 

x 


u 

G 

Ah 


<N 

O 

rs 

Ah 

H 

U 


C  4-1 

o  o 

X  41 
4»  0) 
<n  k- 
G  .c 

^  4-i 

0)  4-> 
J*  01 

g  u 

X 
X  44 
QJ 

U  0) 

c  ^ 

01  '  ‘ 
Cfl 

<u  g 

U  X 
A 

C 

CO  O 

G  *-* 
0  <0 


XJ  X 


xJ  A 
G 

3  XI 
cr  x 
QJ  3 

rH 
0)  4-4 

x: 

H  <0 


C0  M 
fl3  O 

00  (0  4-4 

G 

Cfl  00  O 

XJ 

PS  01 

h  ji  u 

X  XJ  U 
(TJ  G 
X  CO  > 
4J  <  ±J 

C  3 
V  O 

C0  • 

1/1  H  TJ 
4)  X  *tH 
Q>  3 
to  x  cr 

•H  W  rl 
r— I 

C  X 

O  -H  QJ 

^  g  x: 

to  cr  u 


H  CO 
0) 

(T)  U 
X  A 
3  fs  O 

44 

X  -U 
0)  CD  -H 
*H  T3 

X  c  CO 

H  P  TJ 
O  G 

kH 

kn 
3 


CD 

JO 

H  G 

*  X 
CO 

•  4-»  CO 

G  1'  0) 
3  X  U 
O  A  kn 
•H  O  O 
U  k  ^ 
X) 

C  CJ 

O  CD  X 
X  XJ  E 
CO  0)  03 
P  Vh  C 
QJ  O  fs 
a  co  xj 

CO  X 
X  XJ 
XJ 

o 

0J 

>  c 


o 

kH 

QJ 

G 


fs 

X  *h 
CO 

O  Cfl  X 
• — c  Of  a> 

S-S  fi 

CD  kn 
41 
G 


CJ 


CD 

> 

„  O 
QJ  00 
00 


G 
A 
X 
W  CD 

CD 


4->  > 


XJ  QJ 
O  X  * — 1 
6  XI  X 
X  X 
■m  3  G 
O  G  ^3 


co  a>  co 

■H  TJ  -H 
U 

4)  C  C 
4JOO 
O  U  X 
(TJ  CO 

kn  X  _  .  . 

<tj  a  e 

-G  QJ 

o  x  a>  6 

X  CO  *H 
Cl  3  O  00 
G  H  G  (D 
w  44  3  CC 


X 

kH 

G 

U 

G 

O 


G 

O 

o 


u 

G 

A, 


CN 

O 

rs 

& 

u 


co 

■w  41 

a  B 

CD  X 
U  GO 
G  CD 

o  u 
u 

u 

QJ  G 

,g  o 


G  4-1 

o  o 

X  41 
QJ  41 
co  kn 
ctj  x 
X 

(D  XJ 
M  CO 
CTJ  kH 
X 
X  4-H 
CD 
XJ  QJ 
C  JC 
QJ  XJ 
CO 

QJ  G 

G  X 

a 

G 

CO  o 
C  w 
o  <o 

X  X 

w  a. 

G 

•  G  X 
CO  cr  -H 
4)  QJ  G 


CO  w 

G  O 

00  (0  44 

TO 

g  oc  o 


QJ 

e 


Ps  4) 
rj  X 
fH  XJ 

G 

X  « 
-U  < 

C 

41 

CO  » 
CO  X 
41  X 
4 

co  x 


X  CA  X 

C  X 
O  X  41 
W  3  j= 
co  cr  w 
X  X 

A.  X  co 

41 

x  co  a 

*H  U 

G  >\  O 

H  ,C)4H 


X  w 
CO  4J  *H 
H  X 
G  G  CO 
H  G  X 
O  G 
U  G 

•  u  a 
G  3  x 
O  CO  QJ 


eg  qj  iu  w  <d  a> 

B  fC  O  rH  rH 

■w  H  G  B  xi 

-•§•3 


O  O  ^  J3  u 


CO 
•  w 
X  QJ 
G  rH 
O  p. 
rH  o 
O  M 
X 
G 

O  QJ 

*H  W 
CO  CD 
u  u 
ai  o 
a  co 
C0  *H 
■H  X 
X 

o 

QJ  4J 

>  c 

•H  *H 

CO 

O  CO 
«H  QJ 

a.  u 
X  G 
QJ  kn 
QJ 

O  G 

•H  QJ 
W  60 
CO  QJ 
•H  X 
P 

QJ  C 
4->  O 
a  4J 
G  co 

J-J  *H 

g  a. 
x 

a  x 

•H 
QJ  G 

G  rH 


e 

G 

G 

>> 

X 

Q 

kn 

QJ 

G 

rO 

X 

QJ 

G 

kH 

<D 

> 

O 

00 


4J  > 


QJ 

qj  e 

(0  *H 
O  00 

S5 


X 

kH 

G 

U 

G 

O 


c 

o 

u 


XJ 

CO 

QJ 

e 

QJ 

CN 

XJ 

CO 

QJ 

CL 

E 

CO 

kH 

X 

A 

£ 

(/) 

c 

<D 

X 

G 

o 

H 

C 

X 

CD 

•H 

G 

o 

V 

00 

00 

CO  x 

•* 

X 

kH 

CJ 

00 

00 

co  x 

C 

QJ 

G 

X 

G 

c 

QJ 

G 

0 

u 

G 

00  o 

• 

XJ 

co 

CJ 

o 

kH 

G 

00  o 

CJ 

XJ 

X 

QJ 

QJ 

o 

XJ 

C  4h 

o  o 
x  a; 

QJ  QJ 
W  kn 
G  JZ 
JO  u 

QJ  U 
U  V) 

G  U 

•H 
X  4-1 

<v 

4J  QJ 

G  x: 

QJ  V 
CO 

QJ  C 

U  «tH 

a 


V) 

< 


rH  X 


C  X 
O  -H  aj 
G  £ 
co  a*  u 

•rH  »h 
AH  (fl 
QJ 

X  G  U 
X  u 
3  O 
X  44 


a  x 
o  a 
x  o 
a  h  * 
x 
c 

O  QJ  - 
X  4J 
CO  QJ 

k4  kH 
QJ  U 
CL  X 
"  CO  x 
X  X 
X 

o 

QJ  W 

>  a 

X  X 
CO 

O  CO 

X  QJ 

fs 

QJ 

C 


a 

fs 

x 

o 

k4 

CD 

G 

fs 

X 

X 

QJ 

E 

QJ 

> 

O 


XJ 

c 

X 

XJ 

X 

QJ  00 

k< 

CO 

O  co 

QJ 

X 

XJ 

00 

G 

G 

4-1  X 

X 

CO 

QJ  10 

Ah 

O 

CO  x 

e 

CO 

X 

X  X 

X 

X  Eh 

3 

X 

kH 

* 

XJ 

CL 

O 

G 

Q> 

C  C 

CN 

G 

kt 

G 

XJ 

O  0 

o 

• 

3 

X  • 

kH 

A 

u 

w  X 

a*  x  c 

3 

X 

G 

CO  XJ 

QJ 

41 

3  O 

X 

QJ 

W 

X  O 

XJ 

rH  X 

G 

a  0 

A« 

2 

4> 

X  XJ 

QJ 

QJ 

£ 

H 

0 

jC 

o 

X 

X 

U 

X  QJ 

u 

2 

X 

H 

G  0 

X 

X 

x  co 

XJ 

X 

X  X 

X 

G 

■9 

QJ 

J3 

G  O 
x  jr 

Z 

QJ 

o  o 

X 

X> 

XJ 

x  ^ 

G 

O 

CJ 


kH 

G 

Ah 


CM 

O 


Ah 

H 


fs  QJ 
x  X.  X 
X  -U  kH  X 


G  X 
O  O 

X  CD 
0)  QJ 
CO  U 

G  JZ 
X  W 

QJ  W 

k4  C0 

G  k4 
x 
X  X 
QJ 

XJ  QJ 

c  JZ 
QJ  4J 
G 

Q»  c 

W  X 

a 

G 

w  o 

C  4J 
O  CO 
X  X 
4J  Cl 

G 

•  G  X 
w  cr  x 
QJ  QJ  G 


G 

x  co  3 

XJ  <  XJ 

C  p 

01  o 

co  • 

CO  X  x 
QJ  X  -h 

QJ  3 

co  jc  c r 
x  co  x 

X 

G  X 
O  X  QJ 
^  O  JZ 
G  cr  xj 


A  X  CO 
QJ 

X  G  U 

X  H 

3  fs  O 
X  JO  x 

X 

X  XJ 
CO  QJ  X 
X  X 
rC  C  W 
H  G  X 
O  C 
U  G 
•  H  a 
c  g  x 
own 


G  4)  X  XJ  CD  QJ 

BX  O  x  fH 
X  X  G  B  XX 

X  ^  X 

<0  x  x  3  3 

QJ  O  O  ^3  X) 


CO 
-  HJ 

X  QJ 
3  x 

o  a 
o 

a  k« 
x 
c 

O  QJ 
x  XJ 

CO  QJ 

k4  kH 

qj  a 

o.  co 

CO  X 
X  X 
X 

o 

41  XJ 
>  G 

X  X 
CO 

O  co 

X  QJ 

QJ  k4 
QJ 

V  C 

X  QJ 
XJ  00 
G  QJ 
X  X 

kH 

41  G 
XJ  O 

u  w 
G  co 

kH  *H 

g  a 

jC 

U  T3 
X 
QJ  3 

X  X 


6 

G 

G 

fs 

X 

O 

u 

QJ 

G 

>> 

XJ 

X 

CD 

a 

kH 

4J 

> 

O 

00 


G 

o 

x  • 
4-1  > 
O  X 
B 

QJ 

qj  e 
01  x 
O  00 
x:  QJ 
:»  ^ 


<N 

X 

u 

G 

CJ 

c 

o 


X 

XJ 

G 

O 

u 


-  •  i  i  'i-l  I-  illirliilllili'rll'I.iJ  .  'iit  ii.lillUifUWlll..  Iti:  I.',. 


ABSTRACT  CARD 


f 

* 


1 


■.'40.HU  iii . ,,  jh.1ili1Ui1li8ii1a.U4J  1-  u;.  wii..i«-ii  Iiriarf.-.".  ;•  .!SUjH:i::wo m4  .  i..i- .:ir rJiiUJlsm »ini)ilelii liiiftBHamHjlHiimil UMbKHI HMl jUj[[iL8i: 


G 

X 

g 

CO 

re 

c 

o 

T3 

1 

o 

u 

•s_^ 

u 

•H 

QJ 

re 

QJ 

L 

■H 

H 

QJ 

6 

re 

re 

QJ 

CO 

S 

P 

’O 

5 

T5 

B 

u 

w 

* 

re 

3 

G 

QJ 

re 

QJ 

e 

W 

a 

•H 

O 

QJ 

X 

•H 

CL 

3 

P 

• 

a 

T3 

iH 

^3 

re 

G 

re 

03 

33 

O 

CL 

re 

•H 

a 

•H 

00 

QJ 

o 

X 

•o 

•H 

Q 

*o 

vD 

x: 

r-H 

X 

M 

u 

01 

o  u  e 

QJ  o 
03  <0  U 
■H 

©GO) 
?K  <0  X 
H  U 
CO  0 

C  M  4-1 

<0  c  o 


cvs 

rl  H  (/) 

e  CL  03 
•HOW 
H  ^ 

0)  T3  01 

L  W 

a  -a 

H  03 

<|  3  *H 
cr  X 


>>  ©  <0 

H  ,0  C 

QJ  *H 

>  to  6 

*H  *H 

tn  3j  h 

o  o  «) 

iH  «H  H 

p  >  p 
x  o 

<D  H  U 

p  o 


°§:o 

o  TJ  rt  — ^ 


O  T3  (fl 

*H  X) 

03  03 
G  *H  0) 
(0  03  QJ 

a  x> 

X  H  -H 
«J  (fl  > 
G  O 
QJ  CO  U 
X  CL 


X5  I  L 

V  U  C  V 


p 

u 

o 

re  S 

• 

3 

•H 

X 

S 

re  re 

a  o 

X  1 

u 

Li  H 

qj  re 

QJ  Li 

•H 

QJ 

B  re 

re  qj 

a 

P  'O 

5  T5  B 

^  "2 

• 

m  3 

C  «  m 

41  2 

p 

•H  0 

41  £  -rj 

a  3 

• 

a 

-U  H 

x  re  q 

© 

X 

U 

a  re 

*H  O 

00 

QJ 

o  ji 

h 

X  X  _ 

H  B  Q 
D  U  W 
X  U  M 
cO  QJ 
N  O  H 


.G  H  X  4J 

mow 

03  H  *0 

cm  O  P  W 

o  l.  a  > 

wow 

0)  to  u  *H 

•H  i— I 

03  G  W  01 

X  ^ 

tH  *J 
(0  0  03 

G  -U  L-t  «H 
<0  C  O 
•H  M 

>,  M  3 
G  03  C  X 
<0  4J  QJ 


aW  «H 

>  a 


M  09 

(O  *d 

u  * 

«H  (0 


20 


•n 

CJ 

< 

■H 

fH 

re 

00 

Li 

re  re 

*5 

B 

P 

re 

G 

O 

G  *H 

re 

H3 

h-4 

O 

QJ 

•H 

•H 

(0  re 

QJ 

Li 

z 

fH 

L 

re 

cn 

> 

P  5S  T3 

re 

a 

5 

QJ 

T) 

re 

re 

re 

X  fH 

vH 

X 

• 

Li 

re 

re 

J3 

qj  re 

> 

a 

UH 

P  T3 

fr 

QJ 

G 

O 

•H 

•H 

re 

B 

XI 

qj  re 

Li 

0)  G  W  <3J 
^  ^ 
tH  -M 

.(0  0  03 

Q  £  4J  4-t  -H 

w  ©  c  o 

»M  «H  +J 

P  pN  MG 

H  1/3  C  ^3 
C/D  <0  M  W 
C/D  C  0)  S  - 
<2  *H  «H  V)  0C 

J  E  a  u]  c 

CJ  *H  0  QJ  *H 
Z  rH  U  03  01 
3D  0)  *D  03  W 
Li  (0(0 
f"  P  TJ  P 


03  03  QJ  00  > 

£  L  T3  -H  o 

(/)  W  (3  ffl  ^ 

-H  P  3  QJ 
C  03  ©D 

(0  «H  U 
C  0  -H  ^ 

O  0)  L» 

QJ  <0  (0 

£  fH  x>  G 

QJ  «H 
01  P>  (0  £ 

3  *H  *H 

O  03  0)  «H 

•H  O  *G  W 

L  H  -H  O 

re  a  >  a. 

>  X  o 

W  L  Li 


Qi  H  H  ■ 

re 

5kO  M 
XD  ^ 

*>  re  i 

*  QJ  Oh  4 
■'-f  -X 
O  (0  **  I 

^  | 
O  * 
fH  re  I 

0)  c  <*■  • 

T3  «H 
O  X  CL 

X  u  H 


C  X=  z 
O  L  ^ 
•H  *H 
oi  E 
U  CO  • 
QJ  CL 
CL  •  CL 
01  SC 
•H  CO 

Q  re  ^ 

5k 

.  o 


•H 

H 

G 

X 

H 

T3 

tH 

QJ 

u 

G 

x: 

G 

Li 

P 

re 

4h 

• 

<H 

C 

0) 

C 

T3 

rH 

o 

> 

0) 

re 

p 

o 

# 

00 

re 

M 

G 

00 

re 

c 

0) 

C 

u 

T3 

•H 

QJ 

QJ 

o 

M 

X 

•H 

X> 

U 

re 

G 

x: 

u 

3 

c 

O 

QJ 

O0  *H 

<T 

re 

L 

Li 

01 

•H 

X 

•H 

M 

CL 

P 

*H 

H 

> 

«H 

re 

re  c  h 
cl  re 

.  H 
03  re  QJ 
XJ  U  X 


01  C  QJ  QJ 

;k  re  x  x* 


M  O  01 
at 

m  c  o  - 

O  O  XJ 
•H  r 
)  L  03  01 
O  C  *H 

H  re  03 

>  CL  X’ 

re  x  .h  ■ 
x  qj  re 
.  at  g 
x  QJ  re 
x 

0t  M  03 


4-1  _ _ 

sO 


<  3  *H  O 
o-  x  o 
•h  H  C 


xj  i 
qj  h  c 

03  0)  QJ  O0 
E  L  0  *H 
01  Qj  C  01  ' 

*h  p  g  w 

C  Cn  *T3 

re  *h  u 
x  *o  -h  >, 
u  w  L 
a)  >>  re  re 
E  H  ^  c 
QJ  -H 

oi  >  re  b 

G  *H  *H 
O  Ifl  41  H 

rl  O  0  OI 
M  fH  -H  M 
(9  CL  >  (L 

>  X  o 

0)  u  u 

QJ  CL  O 

Uh  Vlh 
^  o  tf)  / 
QJ  (0/ 
*m  G  O  mV 
O  O  T3  C0  N 
•  *H  *0 
3  H  03  O) 

O  C  *H  0) 
•H  re  0)  QJ 

>  a  >v»o 
re  x  h  ^ 
x:  w  re  > 

-  qj  co 
x  qj  re  w 


WV  '  * 

C  QJ 

x 

Li  *  iH 

QJ 

m  q  _c  c 

G 

0) 

X 

Li 

• 

•H 

QJ 

M  C 

x: 

G 

O  > 

re  H 

h  h  re 

O 

> 

*• 

C0 

H 

H 

H 

re 

P  *H 

*  QJ 

P  U  t*-' 

• 

rH 

C  QJ 

B.-H 

«  QJ 

P 

O 

VM 

•  rH 

G  QJ 

(0  01 

y-y  Jxi 

<;  o 

T3 

fH 

o  > 

re 

0) 

C1  rX 

O 

X  fH 

0  > 

•• 

QJ  O 

3?  re 

*2 

QJ 

re 

P  O  •  W3 

QJ 

o 

3d  re 

« 

P 

QJ  <0 

P  o 

• 

00 

3  fH 

CsJ  H  03 

U 

3  00  0)  C 

? 

«H 

^  k! 

CN 

H 

O) 

4J 

3  00 

03 

c 

P 

O  CO  03 

G 

M 

XJ  *H 

P 

O 

CO 

© 

C  M 

T3 

fH 

r-t  X 

rH  (0 

Kcq  41 

QJ 

o 

H  X  *H  03 

tH 

X 

fH  re 

r- 

QJ 

QJ  O 

M  X 

•H 

X3 

re  w 

QJ  G 

<1-  <  o 

0) 

c 

x:  o  3  G 

re 

CjJ 

QJ  c 

< 

(J 

01  G 

X  v 

3 

G 

> 

T3  -h 

o 

QJ 

oo  *h  cr  re 

> 

X)  *H 

O 

QJ 

00  *H 

cr 

re 

re 

O  X 

p  Li 

Li 

CD 

H  X  *H  L» 

re 

o  x: 

P 

Li 

u  re 

fH  X 

•H 

M 

53 

X  u 

H  P 

P  *H 

H  ^  tH  0) 

7Z 

S  o 

H 

P 

P  *H 

iH  5 

re 

m 

USS 

0)  ft 

U  00 


c 

o 

u 

a» 

£ 


c  <m 

o  o 

•V  43 
41  43 
co  k 
cd  js 
.ft  x» 

<u  x> 
k  cn 
03  k 

ft 

73  M 
43 

XJ  43 

ft  £ 


tn 

co 

bO  Cl) 
nj 

03  00 

>v  43 

ft 

ft  XJ 
03 

ft  « 
XJ  < 
C 
43 

to  • 


OJ 

w  ■£ 
ft  co 

ft  73 
O  *H 
XJ  3 

to  cr 

ft  *h 
ft.  ft 


sji 

0  H 


O  • 
xj  73 


xj  73  <o 


ft  >  ft 
*-»  ft  *H 
CO 

co  O 
4)  ft 
O 
k 


>N 

ft 


e-s  e 


CN 

XJ  S 

a  8 

co 

g 

73 

4)  ft 

CO 

0 

k 

u  oo 

bO 

co  kj 

<0 

a  a) 

to 

0  k 

u 

co 

60  O 
XJ 

>,  4) 
H  ft 


3 
T3  O 

k  ft 

ro  u 


(0 

XJ  CO 
0)  4) 
ft  U 
ft.  k 

o  o 

k  vm 


»-t 

w  u 
k  ft 

mJ 

ft 

JD 

w 

VM 

w 

43 

w 

> 

ft 

a. 

vm 

U 

C 

O 

XJ 

c 

73 

X3 

M 

43 

bo 

XJ 

k 

CO  0 

co 

43 

•H 

X* 

bO 

k 

CO 

a  xj 

•h 

73 

co 

43 

co 

cO 

PU 

o  CO 

42 

C 

CO 

•h 

73 

•H 

Pm 

ft  *H 

H 

3 

73 

k 

* 

xj  a 

O 

C 

<D 

C 

c 

CN 

k 

CO 

XJ 

o 

o 

CN 

o 

• 

J  73 

. 

k 

ft. 

U 

XJ 

ft 

• 

O 

r^. 

CO 

cr  -H 

d 

3 

>c 

CO 

CO 

x> 

3* 

r» 

4) 

4)  D 

o 

a) 

43 

k 

ft 

O 

M 

XJ 

H 

ft 

CO 

PX 

e 

(X 

c0 

43  vm 

XJ 

43 

03 

£ 

03 

H 

e 

£ 

o 

ft 

ft 

U 

73 

03 

e 

H 

ft 

H  <0 

8 

£ 

£ 

ft 

CO 

*h 

CJ 

XJ 

£ 

£ 

43 

3 

o 

bO 

CJ 

CO 

VM 

VXJ 

3 

3 

ft 

£ 

43 

pt 

43 

O 

O 

£  £ 

XJ 

VM 

3  « 

2 

73  M 
43 

XJ  0) 

8 -5 
(0 

43  3 
k  *H 

a 

c 

CO  o 
C  t-* 
O  CO 
•H  «H 
XJ  ft. 

03 

•  3  7J 
co  cr  ft 
03  41  3 

xj 

03  4) 

B  £ 
ft  H 
xj 
to 
0) 


i-t  3 
co  cr 
•h  ft 

Cl  r-t 

73  CO 

ft 

3 

•H  £ 
CM 

73 
CO  43 
*H  73 

£  ft 
H  3 
O 
k 
•  k 

e  o 

o  CO 

ft 

XJ  43 
O  ft 
0  ft 
£ 
cm  3 
O  ft 


£  > 
xj  ft 

CO 

CO  O 
43  *H 

o  a 

k  X 
O  4) 


ft 

ft  £ 


CO  73 
43  43 


(0 
CO  ft 
•d  k 

§43 
XJ 
(X  O 
X  <4 
43  k 
«0 

43  £ 
ft  CJ 
ft 

■3.2 

rO  ^ 


43  00 

00 

43  (0 

73  ft 

ft  a 

o  o 

U  ft 

co  k  ; 

■H  g  I 

a  6 


73 

k 

cd 

U 

a 

o 


fu  CJ  M 

1  «rt 

73 

ft 

CO 

3  73 

73 

2  c 

CJ 

XJ 

c 

VM 

XJ 

< 

k  C 

a 

XJ 

: 

3  0  43 

ft 

c 

o 

0 

C 

d  O  43 

d 

0  ft  XJ 

B 

0 

43 

O  ft  XJ 

B 

o 

4 

CD  43 

f0 

o 

73 

43 

CO 

• 

CO  43 

CO 

o 

73  k  k 

u 

v_/ 

43 

43 

03 

t— i 

73  k  k 

P 

VM 

•  2 

H  4|  U 

CO 

k 

oJ 

r— t 

ft  43  o 

3  ft.  CO 

73 

to 

X 

43 

3  ft.  (0 

73 

| 

cr  co  *k 

o 

£ 

XJ 

CO 

cr  co  *h 

O 

ft  ft  73 

k 

ft 

CO 

ft  ft  73 

k 

ft  73 

03 

4> 

XJ 

ft  73 

43 

o 

43 

s 

CO 

c 

-a 

0 

CO 

4)  43  XJ 

{ 

k 

O 

ft 

43  43  X) 

5 

4) 

B 


73  43 
ft  CO 
3  O  00 
ft  ft  43 
CM  3  C* 


k 

03 

&4 


OS 

o 

rx 

<r 

Pm 

H 

u 

§ 


CO 

fs 

CO 

/f. 

£ 

XJ 

43 

E 

43 

CN 

XJ  43 

43 

CN 

"2 

ft. 

0 

CO 

k 

£ 

ft.  6 

CO 

k 

£ 

- 

43 

•H 

<0 

O 

H 

C 

73 

41  ft 

cd 

o 

H 

ft 

73 

O 

00 

oO 

to  vm 

*  *H 

k 

CJ  bO 

&0  co 

VM 

*  ft 

k 

j£ 

c 

43 

<0 

CO 

(0 

a  at 

<0 

CO 

CO 

“ 

o 

k 

cd 

00  0 

• 

XJ  co 

O 

0  k 

co  oo 

O 

• 

XJ  CO 

CJ 

u 

XJ 

73 

43  <13 

■» 

CJ 

XJ 

73 

43  4) 

* 

k 

43 

P 

»-t  a 

a 

k 

43 

ft 

ft  CJ 

d 

■j 

43 

3 

rk 

JZ  73 

O 

&  k 

o 

4)  P 

£ 

73 

O 

ft.  k 

o 

J 

£ 

o 

r-t 

XJ  k 

rH 

O  0 

£  O 

rH  XJ 

k 

ft 

O  O 

X) 

VM 

CO 

CC 

U 

k  Vm 

• 

XJ  «M 

cO 

CO 

u 

k  VM 

• 

•ft 

CO  5 

73 

73 

ft  co 

73 

73 

il 

fi 

VM 

XJ 

<  k 

d 

U 

XJ 

C  *H 

k  < 

k 

ft 

CJ 

XJ 

7 

O 

0 

d 

3 

o 

43  ft 

C 

o  o 

ft 

3 

o 

d 

? 

43 

O 

•H 

k  0 

O 

43 

O 

ft 

^  § 

o 

73 

43 

CO 

CO 

43  3 

CJ 

73  43 

co  • 

CO 

43  CO 

CJ 

> 

43 

43 

CO 

•H  73 

k 

k  d 

VM 

43  43 

CO  ft 

73 

k 

k  ft 

VM 

to 

k 

43 

ft  ft 

43 

u  >, 

CO  k 

43  ft 

4) 

O  >v 

f. 

43  ft 
£  U 

43  XJ 
k  CO 
43  k 
*H 
73  M-t 
4) 

XJ  43 

c  ft 

43 

(0 


43 
X 
ft  CO 

ft  73 

O  ft 

xJ  3 
co  cr 

ft  ft 
p.  ft 

t3  co 


3  Cl 
O'  CO 


■H  «rt  73 


CO  73 
O 


k 
43 

O  CO 
XJ 

a  >% 

-H  ,£3 

CO  73 
03  Q J 


43 

C 

d 

>>  0 

43 

k 

43 

k 

ft 

ft 

^3  «M 

43 

> 

ft. 

VM 

O 

ft 

O 

a 

73  XJ 

ft 

43 

00 

(0 

0 

CO 

43  ft 

XJ 

00 

ft 

XJ 

ft 

73 

CO 

43 

CO 

O 

CO 

£ 

ft  (0 

ft 

73 

ft 

•rt 

ft 

H 

3  73 

k 

XJ 

a 

O  ft 

43 

ft 

ft 

(0 

k  <0 

XJ 

O 

c 

• 

P 

73 

• 

k  a 

U 

x> 

ft 

* 

CO 

V  ft 

d 

3  X 

cfl 

ro 

XJ 

> 

43 

4) 

3 

o 

CD  43 

k 

•H 

o 

w 

XJ 

ft 

d 

ft 

6 

43  Cm 
£  CO 


XJ  43 
O  ft 
0  ft 

t*4  *3 

O  ft 


43  £  43 

H  U  13  41  0 

ft  *H  CO  ft 

.a  03  3  O  00 
3  £  H  £  4) 

O  XJ  4-  J  ftl 


XJ 

k 

to 

P* 


CN 

O 

r- 

(X 

H 


cQ  43 

JO  X3 

4)  xJ 
k  CO 
(0  k 

ft 

*d  •x-i 
43 

X-  03 

C  £ 
43  xj 
CO 

4)  C 
k  -H 

a 

c 

a)  o 
C  4J 

o  co 

ft  ft 

xj  a 
CO 

•  ft  73 
CO  O’  ft 
43  43  3 
XJ  ft 
<0  43  k 

Sg  « 

«J 

CO  4-J 

43  O 


4) 
ft 

■H  CO 
ft  T3 

O  fj 
XJ  P 

co  cr 

ft  *H 
Q.  ft 

73  CO 

ft 

3  ?*> 
ft  ft 

VM 

73 

CO  43 

ft  T> 

6  § 
O 
k 
•  k 

e  ^ 
o  « 

ft 

xj  43 
O  ^ 
S  ^5 

VM  *§ 

O  £ 


o  a 
cr  co 

ft  *h 
•H  73 

43  43 
£  > 
XJ  *H 

co 

to  o 

4)  ft 
9  P? 


(0  73 
ft  o 
73  k 
43 

O  «0 

XJ 

e  >s 
•H  x> 

C0  73 
43  43 


CO 

CO  ft 

73  k 
C  43 

co  xj 

43  k 
*0 

43  £ 
ft  a 
£ 

-D  xt 


43  00 
t>0 

4)  CO 
73  f! 

C  C 

o  o 

XJ  ft  • 
CO  XJ  > 

•MOM 

a  6 

03 

73  43  0 

ft  01  ft 
3  O  00 
ft  £  03 
cm  5  pi 


. .  ''“"f 


ABSTRACT  CARD 


CM  X 

CN 

X 

43 

43 

T3 

43 

a! 

•d 

E 

4> 

La  E 

N 

0) 

■V 
n , 

Li 

N 

01 

cd  5 

M-l 

43 

cd 

o 

*H 

4) 

O  Li 

0) 

AJ 

u 

Li 

0) 

AJ 

W  IM 

cd 

CN 

V— 

<AH 

<d 

AJ 

u 

AJ 

u 

♦  -d 

a 

•H 

TJ 

e 

*H 

•  03 

43 

4) 

4) 

d 

•  a 

M 

□ 

CL 

Li 

c 

•  o 

03 

*rA 

• 

P, 

• 

O 

43 

*H 

• 

•  rH 

Mj 

'V 

Lj 

r— H 

UH 

•  43 

cd 

03 

n 

4) 

Cm 

cd 

43 

•  > 

■H 

AJ 

•H 

tJ 

> 

‘H 

AJ 

•  43 

T3 

cd 

VM 

4) 

T3 

<d 

<A-i 

•  TJ 

*d 

0) 

g 

*d 

d 

01 

O 

*H 

O 

•H 

*  03 

5 

03 

AJ 

43 

3 

01 

AJ 

03  L< 

AJ 

•H 

cd 

0) 

Li 

U 

tH 

cd 

U  T3  5  0) 
O  0)  aj 


U  C  43 
<a*  43  43  E 


01  O  0)  iH 

C  U  U  3 

o  60  cr 

ft  <o  <d  a) 
aj  aj  u 
m  ni  4i 
3  *d  JC  60 

cr  H  c 


e  *h  x:  u 
o  >  lj  nj 

*H  tH 
AJ  *d  >  01 

O  43  '  AJ 

e  u  w  c 

O  G  43 
4)  41  E 
O  -i  £  43 
rH  41  ^ 
03  O  0)  «H 
C  o  L.  3 

o  60  cr 
vi  <q  <q  oi 
aj  aj  u 

CO  <0  03 
3  3  JS  60 

cr  He 

0)  rH  -H 


cn  C  cn  a 

4)  43  41  0) 


0)  C  0)  O 
03  QJ  0)  03 


Li  >  Cd 
03  03  x: 
cl  tj  aj 


V-  >  <d 

a;  03  x: 

a*t3  aj 


*d 

03 

Li  E 

N 

01 

<d  o 

41 

U  Li 

01 

AJ 

vj  U-l 

cd 

AJ 

u 

•  T3 

c 

*tH 

•  43 

43 

'V 

•  a 

Li 

c 

*  0 

43 

•rH 

• 

•  rH 

MJ 

T5 

•  03 

U-i 

co 

41 

•  > 

tH 

AJ 

•H 

•  43 

T? 

cd 

<A-i 

•  T3 

'd 

0) 

O 

•H 

•  4) 

cn 

AJ 

01  Li 

Xj 

<0 

'O  43 

M 

0) 

*H  Jk 

X2 

AJ 

3 

AJ 

03 

•O  Q 
03  e 
0) 

U  VM 

U  03  O 

03  a 

43  1/1  O) 

C  *H  C 
4>  O  O 
O  *H 

>.  4J 
(fl  H  OJ 

C  4)  3 

o  >  or 
a.  *h  4> 
cd  oi 

0)  O  4) 
3  •“»  0) 
a  4) 
*-h  X  x: 
«  u  H 
> 
cd 


■H  X  34 

>  aj  aj 

*H 

13  U  w 

43  AJ 

AJ  AJ  C 
a  C  43 
43  43  E 

H  E  4) 

rH  43  U 
O  4)  *H 
<J  U  3 

60  cr 

<d  oj  q; 
AJ  Li 

<0  43 
TJ  X2  60 
H  C 


AJ  .  (fl 

C  w  u 

03  a 3  0) 

e  o 

^  -H  W 
Li  >  cd 
43  43  x: 
a,  T3  u 


o 


u  o 
4)  B 

to 

Li  UH 

w  43  O 

4)  a 

AJ  CO  01 

C  ^  c 
43  Q  O 

U  *H 
>v  AJ 
01  H  t(J 

C  43  3 
O  >  KT 
CL  -H  43 
TO  01 
43  O  43 
H  (/) 

cl  a) 
^  x  x: 
to  w  h 

cd 

2 


43 

N  0) 

•H  03 
01  AJ 

co 

AJ  O 
C  *H 
4)  *d 
.  Li  C 

4)  *H  • 
U-i  *d 
U-i  (0  43 
*H  AJ  *rl 
*d  CO  >*n 

“O  0) 
O  -H 

Ben  aj 
*h  cd 
x:  cn 
x:  aj 

AJ  43 

■H  £  L 

>  aj  <o 

TJ  5  0) 
43  AJ 
AJ  AJ  C 
U  C  43 
4)  43  E 
•H  E  43 
H  (D  M 
O  03  -H 
U  U  3 
60  cr 
CO  CO  01 
AJ  Li 
Cd  43 

•dX  M 

-  N  e 

r- 1  *»H 

Cd  *H 
aj  •  cd 
C  «  U 
43  43  0) 
6  O 
*H  -H  AJ 

m  >  <d 

03  43  x: 

CL  TJ  AJ 


O 


i  ii:i  i;l|ilf.  i  'i ..  I  ni  l  -UMjiiiHLiMsji  mu  i.siJ  iiii^iiili  >. >l«>nTur>ttiKi[pti|l9U#it^ii{il]f5r>f  AUJfllBU^EIUSIlliiEP 


